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Abstract 

By electron or hole doping quantum antiferromagnets may turn into 
high-temperature superconductors. The low-energy dynamics of anti- 
ferromagnets are governed by their Nambu-Goldstone bosons - - the 
magnons — and are described by an effective field theory analogous 
to chiral perturbation theory for the pions in strong interaction physics. 
In analogy to baryon chiral perturbation theory — the effective theory 
for pions and nucleons — we construct a systematic low-energy effective 
theory for magnons and electrons or holes in an antiferromagnet. The 
effective theory is universal and makes model-independent predictions 
for the entire class of antiferromagnetic cuprates. We present a detailed 
analysis of the symmetries of the Hubbard model and discuss how these 
symmetries manifest themselves in the effective theory. A complete set 
of linearly independent leading contributions to the effective action is 
constructed. The coupling to external electromagnetic fields is also in- 
vestigated. 
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1 Introduction 



Almost 20 years after the discovery of high-temperature superconductivity in layered 
cuprates [1], identifying the dynamical mechanism behind it remains one of the great 
challenges in condensed matter physics. Ordinary low-temperature superconductors 
are weakly coupled electron systems in which phonon exchange mediates an attrac- 
tive interaction that can overcome the Coulomb repulsion between electrons. As 
massless Nambu-Goldstone bosons of the spontaneously broken translation symme- 
try, phonons provide a natural mechanism for Cooper pair formation at low energies 
which is successfully quantified in BCS theory. In contrast to ordinary superconduc- 
tors, layered high-T c cuprates are systems of strongly correlated electrons to which 
the weak coupling BCS theory is not readily applicable. Furthermore, the high 
transition temperatures of cuprate superconductors and the smallness of the isotope 
effect suggest that mechanisms other than phonon exchange may be responsible 
for Cooper pair formation. Since high-temperature superconductors are antiferro- 
magnets before doping, it is natural to suspect (but not generally accepted) that 
magnons — the Nambu-Goldstone bosons of the spontaneously broken SU(2) S spin 
symmetry — may be important for binding electrons or holes into preformed pairs. 

Even if spin fluctuations were not the key to explaining high-temperature super- 
conductivity, the dynamics of charge carriers in an antiferromagnet is an interesting 
topic in itself. There is a vast literature on this subject. The dynamics of holes 
in an antiferromagnet has been investigated, for example, in [2-36]. Understanding 
the dynamics of even just a single hole propagating in an antiferromagnet is a chal- 
lenging problem. One can gain qualitative insight from a picture in which holes hop 
from site to site, leaving a string of flipped spins behind and thus locally destroying 
the antiferromagnetic order. Since the string costs energy proportional to its length, 
one might expect the holes to even be confined and thus have infinite mass. How- 
ever, the locally destroyed antiferromagnetic order may be healed by appropriate 
hole hopping which renders the hole mass finite [7]. Angle resolved photoemission 
spectroscopy experiments [37-40] as well as a number of theoretical investigations 
[6,7,17,34] indicate that the minimum of the dispersion (i.e. of the energy) of a 
single hole corresponds to lattice momenta (±| , ±|) in the Brillouin zone. 

As one adds a second hole, the situation becomes more controversial. For ex- 
ample, there seems to be no consensus on the question if a pair of holes can form a 
bound state or not. If it can, the condensation of such pairs would provide a potential 
mechanism for high-temperature superconductivity. The effective theory to be con- 
structed here can be used to analytically calculate the long-range magnon-mediated 
forces between holes using perturbation theory. It is very interesting to ask what 
happens when one dopes an antiferromagnet with a non-zero density of holes. At 
sufficient doping, experiments show that high-temperature superconductivity may 
arise. It has been argued on theoretical grounds that even an infinitesimal amount 
of doping may affect the antiferromagnetic phase and turn it into a spiral phase [6] . 
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A systematic investigation of this question is also possible using the effective theory 
of this paper, but it will require the use of non-perturbative methods. 

The standard models for antiferromagnets and high-temperature superconduc- 
tors are the Hubbard and t-J model. Since these models are strongly coupled, they 
are not accessible to a systematic analytic treatment. As a consequence, analytic 
calculations in Hubbard-type models usually involve some uncontrolled approxima- 
tions. Unfortunately, due to a severe fermion sign problem, away from half-filling 
these models can currently also not be simulated reliably. Hence, although they 
may indeed contain the relevant physics, Hubbard-type models have not yet led to 
a quantitative understanding of high-T c materials. An alternative to a microscopic 
description using Hubbard-type models is provided by phenomenological models for- 
mulated directly in terms of magnon and electron or hole fields [11,13,25,32,36]. 
Although they may provide qualitative insight, such models do not lead to unam- 
biguous predictions. In this paper, for the first time we introduce a systematic 
low-energy effective field theory for magnons and charge carriers in an antiferro- 
magnet. Based only on symmetries and their spontaneous breakdown, the effective 
theory makes universal predictions for the entire class of antiferromagnetic cuprates. 
Although the effective theory is not renormalizable, it yields unambiguous results 
in a systematic low-energy expansion. In each order of the expansion, the results 
depend only on a finite number of material specific low-energy parameters whose 
values can be determined experimentally. The effective theory is not based on a 
specific microscopic model Hamiltonian but is universally applicable. Furthermore, 
and most important, in contrast to the strongly correlated electrons of Hubbard- 
type models, the electrons and holes of the effective field theory are quasi-particles 
that are weakly coupled to the magnons. Consequently, one may expect that the 
effective theory is more easily solvable than the underlying microscopic models. 

Possible basic applications of the effective theory to be constructed in this paper 
include magnon-magnon, magnon-hole, and magnon-electron scattering as well as 
the determination of long-range magnon-mediated forces between the charge car- 
riers. More ambitious applications could aim at a quantitative explanation of the 
Mott insulator state, the reduction of the staggered magnetization upon doping, the 
formation of a spiral phase, or at a systematic study of potential mechanisms for 
the preformation of electron or hole pairs in the antiferromagnetic phase. When 
such pairs condense they may become the Cooper pairs of high-temperature super- 
conductivity. Except for a derivation of the dispersion relation of charge carriers, 
in this paper we do not consider applications yet, but concentrate entirely on the 
construction of the effective theory itself. 

The construction in this paper is inspired by similar developments in the theory 
of the strong interactions. In contrast to the high-T c problem, where the choice 
of a microscopic model is controversial, there is general agreement that Quantum 
Chromodynamics (QCD) provides the correct microscopic description of the strong 
interactions. Still, similar to Hubbard-type models, solving QCD is notoriously 
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hard. At "half-filling", i.e. in the filled quark Dirac sea that represents the QCD 
vacuum, the SU(2) L <g> SU{2) R chiral symmetry of massless up and down quarks is 
spontaneously broken to the isospin symmetry SU(2)l = r, resulting in three mass- 
less Nambu-Goldstone pions. This is analogous to the spontaneous breaking of 
the SU{2) S spin symmetry down to U(l) s that leads to antiferromagnetism. The 
corresponding Nambu-Goldstone bosons — in this case two magnons — are thus 
analogous to the pions of the strong interactions. It is possible to study chiral sym- 
metry breaking in the QCD vacuum in numerical simulations of lattice QCD, just as 
it is possible to study antiferromagnetism by simulating the Hubbard model at half- 
filling. However, it is very useful to also investigate these phenomena with effective 
field theories. The low-energy effective theory for pions was pioneered by Weinberg 
[41] and formulated as a systematic expansion in Gasser's and Leutwyler's chiral 
perturbation theory [42]. Based on symmetry considerations and the observation 
that chiral symmetry is spontaneously broken, chiral perturbation theory makes rig- 
orous predictions about the pion dynamics in terms of a few low-energy parameters 
such as the pion decay constant, the chiral condensate, and the Gasser-Leutwyler co- 
efficients. Once these parameters are determined, either experimentally or through 
lattice QCD calculations, the effective theory makes unambiguous predictions in the 
low-energy domain. 

Chiral perturbation theory can be applied to any Nambu-Goldstone phenome- 
non, and has indeed been used for both ferro- [43, 44] and antiferromagnetic magnons 
[45-50]. To lowest order, for antiferromagnetic magnons the low-energy parameters 
of chiral perturbation theory are the spin stiffness p s and the spin-wave velocity 
c. At low energies chiral perturbation theory describes all aspects of the magnon 
dynamics just in terms of these two parameters. For example, the low-energy physics 
of the Hubbard model at half-filling is completely described by the effective theory 
once p s and c have been determined in terms of the Hubbard model parameters t 
and U . 

A numerical challenge in high-T c physics is to simulate the Hubbard model away 
from half-filling. This requires a solution of the corresponding fermion sign problem. 
Similarly, simulating lattice QCD at non-zero baryon chemical potential, i.e. after 
"doping" the QCD vacuum with quarks, is prevented by a severe complex action 
problem. Like for high-T c materials at sufficient doping, one expects that QCD at 
sufficiently high baryon density becomes a superconductor, in that case for the color 
charge carried by quarks and gluons [51]. In contrast to high-temperature supercon- 
ductivity, the mechanism responsible for color-superconductivity is well understood 
in terms of one-gluon exchange. Color-superconductivity requires very large baryon 
densities and may thus arise only in the core of compact neutron or quark stars. 
However, superconductivity — not of color but of ordinary electric charge — is also 
known to exist at more moderate baryon densities. In particular, pairing of protons 
or neutrons inside large nuclei or neutron stars leads to superconductivity or su- 
perfluidity. Understanding the mechanism of nucleon pairing from the microscopic 
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QCD theory may be as hard as understanding the mechanism for high-temperature 
superconductivity directly from the Hubbard model. Instead it is much more useful 
to employ a systematic low-energy effective theory whose parameters can be deter- 
mined from the underlying microscopic physics. In nuclear physics effective field 
theory has recently led to some progress in describing the forces between nucle- 
ons in terms of just a few low-energy parameters [52-61], while phenomenological 
models involve a much larger number of adjustable parameters. Also steps towards 
describing nuclear matter with effective field theories have already been taken [62- 
66] . The goal of the present paper is to develop a similar effective theory describing 
the interactions between the charge carriers in an antiferromagnet through magnon 
exchange. Remarkably, some physical phenomena that are practically inaccessible 
to microscopic Hubbard-type models even by numerical simulation can be tackled 
analytically in the effective field theory framework. 

An ambitious goal of the effective theory approach is to systematically investigate 
possible mechanisms for the preformation of electron or hole pairs as a potential step 
towards understanding high-temperature superconductivity. It is an experimental 
fact that antiferromagnetism is destroyed before one enters the superconducting 
phase. How can magnon exchange then possibly provide a mechanism relevant for 
Cooper pair preformation? The destruction of antiferromagnetism just means the 
absence of infinite-range antiferromagnetic order. Antiferromagnetic correlations, 
although only of finite range, exist even in the superconducting phase. The finite 
correlation length implies that the magnons have developed a massgap, but they 
may still exist as relevant low-energy degrees of freedom. In particular, in 2 + 1 
dimensions, as a consequence of the Hohenberg-Mermin-Wagner-Coleman theorem, 
magnons pick up a mass that is exponentially small in the inverse temperature 
[67,68]. The generation of the massgap is a non-perturbative phenomenon that is 
well within the applicability range of the effective theory, although infinite-range 
antiferromagnetic order exists only at zero temperature. Similarly, an effective the- 
ory for magnons and electrons or holes remains valid in the superconducting phase 
as long as the magnons remain among the lightest degrees of freedom. Again, this 
is similar to QCD where pions are not exactly massless either — in that case as a 
result of explicit chiral symmetry breaking due to non-zero quark masses. Although 
pions are hence only pseudo-Nambu-Goldstone bosons, chiral perturbation theory 
remains perfectly well applicable. 

The low-energy effective theory for magnons and charge carriers to be developed 
here is the condensed matter analog of baryon chiral perturbation theory in strong 
interaction physics [69-73]. The effective theory is based on a non-linear realization 
of the spontaneously broken symmetry [74, 75]. The terms in the low-energy effective 
Lagrangian are organized according to the number of derivatives they contain. The 
lowest energy physics is dominated by the terms with the smallest number of deriva- 
tives, while effects at higher energies are taken into account systematically through 
higher- derivative terms. A key ingredient in constructing the effective Lagrangian 
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are symmetry considerations. At a given order of the low-energy expansion, i.e. for 
a given number of derivatives, all terms consistent with the symmetries must be 
included in the effective Lagrangian, with a low-energy parameter that determines 
the strength of the corresponding interaction. For cuprates the most important 
symmetries are the SU(2) S spin symmetry which is spontaneously broken down to 
U(l) s in the antiferromagnetic phase, as well as the U(1)q fermion number symme- 
try whose breakdown signals superconductivity. Other relevant symmetries include 
translation by one lattice spacing which changes the sign of the staggered magneti- 
zation, 90 degrees rotations and reflections of the square crystal lattice, as well as 
time-reversal. In addition to these generic symmetries of high-T c materials, the Hub- 
bard model possesses an SU(2)q symmetry discussed by Yang and Zhang [76,77] 
which is a non-Abelian extension of the charge symmetry U(1)q. This symmetry is 
not expected to be present in generic cuprate materials, but may still be a relevant 
approximate symmetry in specific samples. 

In this paper we ignore phonons, assuming that they do not play an important 
role for high-temperature superconductivity. For example, in the Hubbard model a 
rigid lattice which does not have its own physical degrees of freedom is put by hand. 
Of course, in the actual high-T c materials a crystal lattice arises as a result of the 
spontaneous breakdown of translation and Galilean (or more precisely Poincare) 
invariance. The corresponding Nambu-Goldstone bosons are the phonons. The 
role of phonons and their possible interplay with magnons can also be investigated 
systematically in the framework of low-energy effective field theory. 

We also consider the coupling of antiferromagnets to external electromagnetic 
fields which can be used to probe the dynamics of magnons and electrons or holes. 
As first noted by Frohlich and Studer, in non-relativistic condensed matter external 
electromagnetic fields E(x) and B(x) enter the dynamics in the form of non-Abelian 
vector potentials for the SU(2) S spin symmetry [78]. We use this observation to cou- 

— * 

pie both the microscopic Hubbard model and the effective theory to external E{x) 

— * 

and B(x) fields. As discussed in detail in [79], the electromagnetic couplings are 
the condensed matter analog of the weak interactions in particle physics. These are 
described by an SU(2) l <S>U(1)y gauge theory, which turns part of QCD's global chi- 
ral symmetry into a gauge symmetry. Remarkably, the electromagnetic couplings of 
non-relativistic condensed matter are described by a local SU (2) s <S>U(1)q symmetry 
which is the condensed matter analog of the SU(2) L <g> U(l)y symmetry in particle 
physics. Some correspondences between QCD and antiferromagnets are summarized 
in table 1. Connections between QCD and condensed matter physics have also been 
discussed in [80]. 

The rest of this paper is organized as follows. Section 2 contains a symmetry 
analysis of the Hubbard model as a concrete example for an underlying microscopic 
system. In section 3 the effective theory for magnons is reviewed and the non-linear 
realization of the SU(2) S spin symmetry is constructed. In section 4 the Hubbard 
model is coupled to a magnon background field. In this way the fields of the effective 
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QCD 


Antiferromagnets 


broken phase 


hadronic vacuum 


antiferromagnetic phase 


global symmetry 


chiral symmetry 


spin rotations 


symmetry group G 


SU(2) L ®SU(2) R 


SU(2) S 


unbroken subgroup H 


SU(2) L=R 


U(l) s 


Goldstone boson 


pion 


magnon 


Goldstone field in G/H 


U(x) G SU(2) 


e(x) G S 2 


order parameter 


chiral condensate 


staggered magnetization 


coupling strength 


pion decay constant F„ 


spin stiffness p s 


propagation speed 


velocity of light 


spin-wave velocity 


conserved charge 


baryon number U(1)b 


electric charge U(l)n 


charged particle 


nucleon or antinucleon 


electron or hole 


long-range force 

o o 


pion exchange 


magnon exchange 


weak probes 


electroweak fields 


electromagnetic fields 


local symmetry of 


electroweak 


local 


weak probes 


SU{2) L ®U{l) Y 


SU(2) a ®U(l) Q 


dense phase 


nuclear or quark matter 


high-T c superconductor 


microscopic description 


lattice QCD 


Hubbard-type models 


effective description 


chiral perturbation 


magnon effective 


of Goldstone bosons 


theory 


theory 


effective description 


baryon chiral 


effective theory 


of charged fields 


perturbation theory 


presented here 



Table 1: Correspondences between QCD and antiferromagnets. 

theory inherit their transformation properties under the various symmetries from 
the underlying microscopic degrees of freedom. In section 5 the low-energy effective 
theory for magnons and charge carriers is developed and the leading terms in a 
systematic low-energy expansion of the effective action are constructed. This section 
also contains an application of the effective theory to the dispersion relations of 
charge carriers. Section 6 treats the t-J model and its effective theory as a special 
case of systems with holes as the only charge carriers. In section 7 the Hubbard 
model as well as its effective theory are coupled to external electromagnetic fields. 
Finally, section 8 contains our conclusions, while some technical details are discussed 
in two appendices. 



2 Symmetries of the Hubbard Model 

In order to have a concrete microscopic system for which we will then construct a 
low-energy effective theory, we consider the Hubbard model. The Hubbard model 
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just serves as one representative of a large class of systems, including the actual 
high-T c materials. Here it is essential that the Hubbard model shares important 
symmetries, e.g. an SU(2) S spin symmetry and a U(1)q fermion number symmetry 
with these materials. In the Hubbard model at half-filling the U(1)q symmetry 
even extends to an SU(2)q symmetry. The SU(2)q symmetry is not exact in actual 
materials, but may still be approximately realized and will also be investigated in 
the framework of the effective theory. 



2.1 Hamiltonian and Generic Continuous Symmetries 

The Hubbard model is defined by the Hamiltonian 

H = ~t £(4t c *-ht + c[^c x} + c\ Mi + cl +h c xl ) 

x,i 

X X 

Here x denotes the sites of a 2-dimensional square lattice and i is a vector of length 
a (where a is the lattice spacing) pointing in the i-direction. Furthermore, t is the 
nearest-neighbor hopping parameter, while U > is the strength of the screened on- 
site Coulomb repulsion, and fx' is the chemical potential for fermion number. The 
fermion creation and annihilation operators obey the standard anticommutation 
relations 

{4a> C Vs'} = Sxy$ss>, {c X s, C ys >} = {c\ s , c\ gl } = 0. (2.2) 

We also introduce the SU(2) S Pauli spinor 

in terms of which (up to an irrelevant constant) the Hamiltonian takes the manifestly 
S'L r (2) s -invariant form 

H = ~t E(4c x+i + cl +fx ) + I Y,(cic x - l) 2 - /iE(4c, - 1). (2.4) 

x,i x x 

Here // = // — \U is the chemical potential for the fermion number relative to half- 
filling, i.e. fx = implies an average density of one fermion per lattice site. The 
corresponding U(1)q symmetry is generated by the charge operator 

Q = E^ = E(4c,-l). (2.5) 

X X 

Again, we count fermion number relative to half-filling. The SU(2) S symmetry is 
generated by the total spin 

£=E& = E4 f c *> ( 2 - 6 ) 
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where a are the Pauli matrices. It is easy to see that the above Hamiltonian conserves 
both fermion number and spin, i.e. [H, Q] = [if, S] = 0, and that [Q, S] = 0. 
The infinitesimal generators S of SU(2) S (which obey the standard commutation 
relations [S a , Sb] = ie a bcS c ) can be used to construct a unitary operator 

V = exp(iff- S), (2.7) 

which implements the corresponding symmetry transformations in the Hilbert space 
of the theory. In particular, the transformed annihilation operators take the form 

— * — * 

c' x = V^c x V = exp(iff- -)c x = gc x , g = exp(iff ■ — ) e SU(2) S . (2.8) 

Similarly, the U(1)q transformations are implemented by a unitary operator 

W = exp(iu;Q), (2.9) 

such that 

Q c x = W^c x W = exp(icu)c x , exp(itu) e U(1) Q . (2.10) 

For large positive U, at half-filling, the repulsive Hubbard model reduces to the 
antiferromagnetic spin | quantum Heisenberg model with the Hamiltonian 

H = jJ2S x -S x+h (2.11) 

x,i 

where the exchange coupling is given by J = 2t 2 /U. This follows to second order 
of perturbation theory int/U. To leading order, i.e. completely ignoring the kinetic 
term proportional to t, there is an enormous number of degenerate ground states. 
Irrespective of spin, any state with exactly one fermion occupying each lattice site 
avoids the on-site Coulomb repulsion and thus represents a ground state for t — 0. 
There is no correction at order t/U. In second order of degenerate perturbation 
theory, a spin can virtually hop to a neighboring site occupied by a fermion with 
opposite spin and then hop back. On the other hand, virtual hops to sites occupied 
by a fermion with the same spin orientation are forbidden by the Pauli principle. 
This favors antiparallel spins and leads to the antiferromagnetic Heisenberg model 
of eq.(2.11). 



2.2 Discrete Symmetries 

Since the Hubbard model at half-filling leads to antiferromagnetism, another impor- 
tant symmetry is translation by one lattice spacing (in the i-direction), which flips 
the sign of the staggered magnetization vector 

m s = ei- 1 )^- ( 2 - 12 ) 

X 
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The factor (— l) x = (— \^ x ^+ x ^)/ a distinguishes between the sites of the even and odd 
sublattice. The points on the even sublattice A have (— l) x = 1 while the points on 
the odd sublattice B have (— l) x = — 1. The displacement symmetry is generated 
by a unitary operator D which acts as 

D c x = D^c x D = c x+h (2.13) 

and for which [H, D] = 0. Obviously, both the U(1)q and the SU(2) S symmetry 
commute with the displacement, i.e. [Q, D] = [S, D] = 0. In the effective theory it 
will be useful to also consider a related symmetry D' which combines D with the 
spin rotation g = ia 2 . This symmetry acts as 

D 'c x = D*c x D' = (ia 2 ) D c x = (ia 2 )c x+l (2.14) 

Also note that [H, D'\ = [D, D'\ = [Q, D'] = 0, but [S, D'] ^ 0. 

In non-relativistic physics orbital angular momentum and spin are separately 
conserved and spin plays the role of an internal quantum number. Indeed, in the 
Hubbard model the SU(2) S spin symmetry is completely independent of the 90 
degrees rotation invariance of the spatial lattice. The 90 degrees rotation O acts on 
a spatial point x = (xi,x 2 ) as Ox = (—x 2 ,xi). Under the symmetry O the fermion 
operators transform as 

°c x = ] c x O = c 0x . (2.15) 

Parity turns x into (— x±, — x 2 ) and is equivalent to a 180 degrees rotation in two 
dimensions. Hence, it is more useful to consider the spatial reflection R at the X\- 
axis which turns x into Rx = (xi,—x 2 ). Under this transformation the fermion 
operators transform as 

R c x = rfc x R = c Rx . (2.16) 

The reflection at the orthogonal X2-axis is a combination of the reflection R and 
the rotation O. One can also consider the reflection at an axis half between lattice 
points. This transformation is a combination of R with the displacement symmetry 
D. Similarly, a reflection at a lattice diagonal is a combination of R and O. An- 
other important symmetry is time-reversal which is implemented by an antiunitary 
operator T. 

It should be pointed out that, unlike the actual high-T c materials, the Hubbard 
model is not Galilean invariant: in the actual materials translation as well as Galilean 
invariance are spontaneously broken by the formation of the crystal lattice. The 
corresponding Nambu-Goldstone bosons are the phonons which are known to play 
a central role in ordinary low-T c superconductivity. In the Hubbard model, on 
the other hand, the lattice is imposed by hand, and thus translation and Galilean 
invariance are explicitly broken. In particular, phonons cannot arise because the 
lattice does not have its own physical degrees of freedom. 
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2.3 SU(2) Q Symmetry 



As first noted by Yang and Zhang [76, 77], at half-filling (i.e. for /i — 0) the Hubbard 
model possesses a non-Abelian extension SU(2)q of the fermion number symmetry 
U(1)q generated by 

Q + = E(-i) x 4t4|, Q- = £(-1)*^, Q 3 = £ i(i TQeT + 4 iC;4 - 1) = 1q. 

X X X Z Z 

(2.17) 

Writing Q ± = Q 1 ± iQ 2 , it is straightforward to show that, for /i — 0, indeed 
[if, Q] = 0. Also the SU(2)q symmetry commutes with the SU(2) S symmetry, i.e. 
[Q a , S b ] = 0, but it does not commute with the displacement symmetry because 
D^Q ± D = -Q ± . For the same reason [Q, D'] ^ 0. 



Introducing the SU (2)q spinor 



dx = ( ) , (2.18) 



which obeys the standard anticommutation relations 

{4a, = S xySab, {4a, ^fe} = {4a> 4^ = ' ( 2 ' 19 ) 

one writes 

— * 

<3 = £<3* = £4? 4, (2.20) 

where a are again the Pauli matrices now operating in SU(2)q space. The infinites- 
imal generators Q of SU{2)q can be used to construct a unitary operator 

W = exp(iw-$), (2.21) 

which implements the corresponding symmetry transformations in the Hilbert space 
of the theory. The transformed SU (2)q spinors are then given by 

— * — * 

Qd x = W ] d x W = exp(itu ■ ^)d x = Qd X: SI = exp(iCu ■ |) e SU{2) Q . (2.22) 
In terms of the SU (2)q spinors the Hamiltonian takes the form 

H = -t£(44+* + d[ +1 d x ) - ^E(44 - 1) 2 -^£4Mz- (2.23) 

x,i 'x x 

The first two terms on the right-hand side are manifestly SU( 2) q- invariant, while 
away from half-filling (i.e. for /i ^ 0) the chemical potential term explicitly breaks 
the SU{2)q symmetry down to U(1)q. 
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Finally, we introduce a matrix-valued fermion operator 



V c -i -(- 1 ) 4 T / 

which displays both the SU(2) S and the SU(2)q symmetries in a compact form. 
The first column of C x is the SU(2) S spinor c x , while the second column is another 
SU(2) S spinor which transforms exactly like c x . The first row of C x is the SU(2)q 
spinor d x , while the second row is another SU(2)q spinor which transforms exactly 
like d x . Under combined SU(2) S and SU(2)q transformations C x transforms as 

$C' X = gC x Q T . (2.25) 

Since the SU(2) S symmetry acts on the left while the SU(2)q symmetry acts on the 
right, it is now manifest that the two symmetry operations commute. Under the 
displacement symmetry one obtains 

D C X = C x+i a 3 . (2.26) 

The appearance of a 3 on the right is due to the factor (— l) x and confirms that the 
displacement symmetry commutes with all SU(2) S transformations, but only with 
the Abelian U(1)q (and not with all SU(2)q) transformations. Similarly, under the 
symmetry D' one finds 

°'C X = (ia 2 )C x+ -a 3 . (2.27) 

The Hamiltonian can now be expressed in a manifestly SU(2) S -, U(1)q-, D-, and 
D'-invariant form 

# = - 1 £ Tr [Ct + C[ f x ] + ^ £ Tr[CtC x CtC x ] - ^ Tr[CtC x a 3 ] . (2.28) 

Z X,i lZ X Z X 

The chemical potential term is only U(1)q invariant, while the other two terms are 
manifestly «S?7(2)Q-invariant. 



3 Effective Theory for Magnons 

Before doping, the high-T c materials are quantum antiferromagnets in which the 
SU(2) S spin symmetry is spontaneously broken down to U(l) s . The low-energy 
physics of antiferromagnets is dominated by the corresponding Nambu-Goldstone 
bosons — the magnons. Chiral perturbation theory, which was originally developed 
for the Nambu-Goldstone pions of QCD, is a systematic low-energy expansion that 
has also been applied to magnons [43-50]. In this section we review the basic 
features of magnon chiral perturbation theory. As a necessary prerequisite for the 
coupling of magnons to charge carriers, we also construct the non-linear realization 
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of a spontaneously broken SU(2) S symmetry, which then appears as a local U(l) s 
symmetry in the unbroken subgroup. This is analogous to baryon chiral perturbation 
theory in which the spontaneously broken SU(2)l ® SU(2)r chiral symmetry of 
QCD is implemented on the nucleon fields as a local SU(2)l = r transformation in 
the unbroken isospin subgroup. 



3.1 Continuous Symmetries of the Effective Action 

The undoped precursors of high-temperature layered cuprate superconductors are 
quantum antiferromagnets. At half-filling, also the Hubbard model displays anti- 
ferromagnetism. In these systems, at least at zero temperature, the spin rotational 
symmetry G = SU(2) S is spontaneously broken down to the subgroup H = U(l) s 
by the formation of a staggered magnetization. The U(1)q symmetry, on the other 
hand, remains unbroken until one reaches the superconducting phase. In the Hub- 
bard model even the SU(2)q symmetry remains unbroken at half-filling but is ex- 
plicitly broken down to U(1)q for /i ^ 0. As a consequence of Goldstone's theorem, 
there are two massless bosons — the antiferromagnetic spin-waves or magnons, 
which are described by a unit- vector field 

e(x) = (e 1 (x),e 2 (x),e 3 (x)), e(x) 2 = 1 (3.1) 

in the coset space G/H = SU(2) S /U(1) S = S 2 . Here t) denotes a 

point in Euclidean space-time. The vector e(x) describes the direction of the local 
staggered magnetization. The leading order terms in the Euclidean action of the 
low-energy effective theory for the magnons take the form [67, 68] 

S[e\ = J d 2 x dt ^ (die- <9,e + ^d t e- <9 t e) . (3.2) 

The index i e {1,2} labels the two spatial directions, while the index t refers 
to the Euclidean time-direction. The parameter p s is the spin stiffness and c is 
the spin-wave velocity. For the antiferromagnetic Heisenberg model of eq.(2.11) 
these low-energy parameters have been determined very precisely in Monte Carlo 
calculations [81,82] resulting in p s = 0.186(4) J, c = 1.68(1) Ja, where J is the 
exchange coupling of the Heisenberg model and a is the lattice spacing. The leading 
terms in the magnon effective action are "Poincare" -invariant with the spin- wave 
velocity c playing the role of the velocity of light. Consequently, antiferromagnetic 
magnons have a "relativistic" spectrum. The "Poincare" symmetry emerges only at 
low energies as a consequence of the discrete lattice rotation invariance. However, 
higher-derivative terms relevant at higher energies are in general not invariant. 

In the following we prefer to work with an alternative representation of the 
magnon field using 2x2 Hermitean projection matrices P(x) that obey 

P( X y = P{x), TrP(x) = 1, P(x) 2 = P(x), (3.3) 
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and are given by 



m = i (1 + e» • 3) = i ( ^ ) . (3,) 



In the above CP(1) language, the lowest-order effective action of eq.(3.2) takes the 
form 

S[P] = J d 2 x dt p s Ti f ~ 



d t Pd t P + -d t Pd t P 



c 2 



(3.5) 



This action is invariant under the global transformations g G SU(2) S of eq.(2.8), 

P(x)' = gP(x)gl (3.6) 

Note that the magnon field P(x) is invariant under the charge symmetries U(1)q 
and SU(2) Q , i.e. $P{x) = P(x). 

3.2 Discrete Symmetries of Magnon Fields 

Under the displacement D by one lattice spacing the staggered magnetization chan- 
ges sign, i.e. 

D e(x) = -e(x) => D P{x) = 1 - P{x). (3.7) 

Let us again combine D with the spin rotation g = ia2, which results in the trans- 
formation D' with 

D 'P(x) = (ia 2 ) D P(x)(ta 2 y = (ia 2 )[l - P{x)\{ia 2 )^ = P(x)*, (3.8) 

reminiscent of charge conjugation in particle physics. 

The Hubbard model is invariant under translations by an integer multiple of the 
lattice spacing. As we have seen, due to the antiferromagnetic order, the displace- 
ment D by one lattice spacing (which connects the two sublattices A and B) plays 
a special role. In particular, in the effective theory it manifests itself as an internal 
symmetry that changes the sign of e(x). Translations by an even number of lattice 
spacings (which do not mix the sublattices), on the other hand, manifest themselves 
as ordinary translations in the effective theory. It should be noted that in the ef- 
fective theory one need not distinguish between the displacement symmetries D for 
the two spatial directions, since they are related by an ordinary translation by two 
lattice spacings (one in the 1- and one in the 2-direction). 

When we decompose a space-time vector x = (x\,x 2 ,t) into its spatial and 
temporal components, the 90 degrees rotation O acts on x as Ox = (—x 2 ,x\,t). 
Under the symmetry O the magnon field transforms as 



o 



P(x) = P{Ox). (3.9) 
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Similarly, under the spatial reflection R at the xi-axis, which turns x into Rx = 
(xi, —x 2 ,t), the magnon field transforms as 

R P(x)=P(Rx). (3.10) 

Had we not treated spin as an internal quantum number, it would also be directly 
affected by the spatial reflection. Since spin is a form of angular momentum, it 
transforms like the orbital angular momentum L = f x p of a particle, which is a 
pseudo-vector and thus changes into R L = (— L ± ,L 2 ,— L 3 ) under the reflection R. 
This is equivalent to a 180 degrees SU (2) s rotation around the 2-direction. Since we 
treat SU(2) S as an exact internal symmetry, the pure spatial inversion R (without 
180 degrees rotation of the spin) is also a symmetry. 

Another important symmetry is time- reversal T which turns x = (xi,x 2 ,t) into 
Tx = (xi, x 2 , — t). In a Hamiltonian description time- reversal is represented by 
an antiunitary operator. Here we discuss time-reversal in the framework of the Eu- 
clidean path integral. Again, the spin transforms like the orbital angular momentum 
L of a particle. The momentum p changes sign under time-reversal and so does L, 
i.e. T L = —L. 1 Consequently, under T the staggered magnetization vector (which 
is built from microscopic spins) transforms as 

T e(x) = -e(Tx) =>■ T P(x) = 1 - P(Tx) = D P{Tx). (3.11) 

Hence, time-reversal is closely related to the displacement symmetry of eq.(3.7). 
Just like the displacement symmetry D, time-reversal is spontaneously broken in 
an antiferromagnet. However, in contrast to a ferro magnet, the combination TD 
of time-reversal and the displacement symmetry remains unbroken. Previously we 
have combined the displacement symmetry D with the SU{2) S spin rotation ia 2 in 
order to obtain the unbroken symmetry D' . In order to obtain an unbroken variant 
T' of time-reversal we now combine T with the spin rotation ia 2 which yields 

T 'P(x) = (ia 2 ) T P(x)(ta 2 y = (ia 2 ) D P{Tx){ia 2 y = D 'P(Tx). (3.12) 

3.3 Non-Linear Realization of the SU(2) S Symmetry 

In order to couple electron or hole fields to the magnons one must construct a 
non-linear realization of the spontaneously broken SU(2) S symmetry which then 
manifests itself as a local symmetry in the unbroken U(l) s subgroup of SU(2) S . 
This local transformation is constructed from the global transformation g G SU(2) S 
as well as from the local magnon field P(x) as follows: one first diagonalizes the 

^^Note that —L does not obey the angular momentum commutation relations. This is a conse- 
quence of the antiunitary nature of T which does not represent an ordinary symmetry (implemented 
by a unitary transformation) in Hilbert space. 
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magnon field by a unitary transformation u(x) G SU(2) S , i.e. 

u(x)P(x)u(xy = ^(l + a 3 )=^ 1 Q jQ, «n(x)>0. (3.13) 

Note that, due to its projector properties, P(x) has eigenvalues and 1. In order 
to make u(x) uniquely defined, we demand that the element u n (x) is real and non- 
negative. Otherwise the diagonalizing matrix u(x) would be defined only up to a 
U(l) s phase. Using eq.(3.4) and spherical coordinates for e(x), i.e. 

e(x) = (sm0(x) cosip(x), sin9(x) smip(x), cos9(x)), (3-14) 

one obtains 

= 1 ( l + e 3 (a;) e 1 (x) - ie 2 (x) \ 

y/2(l + e 3 (x)) V -ei(x) - ie 2 (x) 1 + e 3 (x) J 

( cos(|0(a;)) sin(|0(x)) exp(— iip(x)) 

\ — sin(|^(x)) exp(itp(x)) cos(|^(x)) 

Under a global SU(2) S transformation g the diagonalizing field u(x) transforms as 

u(x)' = h(x)u(x)g\ u u {x)'>0, (3.16) 

which implicitly defines the non-linear symmetry transformation 

m = = ( eXP(l ; (l)) exp( _° ia(l)) ) 6 U(l),. (3.17) 

The transformation /i(x) is uniquely defined since we demand that Un(x)' is again 
real and non-negative. Note that with this definition of h(x) indeed 

u(x)'P(xy U (xf = ±(l + a 3 ). (3.18) 

Interestingly, the global SU(2) S transformation g manifests itself in the form of a 
local transformation h(x) G U(l) s which inherits its x-dependence from the magnon 
field P(x). 

We still need to show that the SU(2) S group structure g = g 2 g\ is inherited by 
the non-linear U(l) s realization, i.e. h(x) = h 2 (x)hi(x). First, we perform the global 
SU(2) S transformation g±, i.e. 

P{x)' = gi P(x)gl u{x)' = h l {x)u{x)gl (3.19) 

which defines the non-linear realization h\{x). Then we perform the subsequent 
global transformation g 2 which defines the non-linear realization h 2 (x), i.e. 

P(x)" = g 2 P{x)'g\ = g 2gi P(x)(g 2 g^ = gP{x)g\ 

u(x)" = h 2 {x)u{x)'g\ = h 2 (x)h 1 (x)u(x)(g 2 g 1 ) t = h(x)u(x)g ] . (3.20) 
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This indeed implies the correct group structure h(x) = h 2 (x)h\(x). 

Under the displacement symmetry D the sign-change of the staggered magneti- 
zation e(x) implies 

D = 1 ( l-e 3 {x) -ei(x) +ie 2 (x) \ 

y/2(l-e 3 (x)) V e 1 (x)+ie 2 {x) 1 - e 3 (x) J 

( sin(|0(x)) — cos(|0(z)) exp(—iip(x)) \ 

y cos(|0(z)) exp(i(p(x)) sin(|0(z)) y 

= r(x)u(a:), (3.21) 



where 



r(x) 



1/^0 —ei(x) + ie 2 (x) 





v /e 1 (x) 2 + e 2 (a;) 2 V ^ + ie *( x ) 

( -exp(-i^(a;)) \ 
I exp(i<p(x)) i ' 



(3.22) 



Note that d t(x) = —t{x) = t(x) J[ , such that 

DD u(x) = D r(x) D u(x) = t(x) ] t(x)u(x) = u(x), (3.23) 

as one would expect for the displacement symmetry. It should also be noted that 
- like the SU(2) S symmetry — the displacement symmetry is also spontaneously 
broken and hence realized in a non-linear (i.e. magnon-field-dependent) manner. 
Similarly, under the displacement symmetry D' one finds D u(x) = h(x) D u(x)g^ 
with g = ia 2 . For this particular g the local transformation takes the form h(x) = 
{i(t 2 )t{x)\ such that 

D 'u(x) = u(x)*. (3.24) 

In contrast to the displacement symmetry D, the symmetry D' is not spontaneously 
broken and is thus realized in a linear (i.e. magnon- field-independent) manner. 

In the next step we consider the anti-Hermitean field 

Vp(x) = u(x)d^u{x)\ (3.25) 
which transforms under SU(2) S as 

v^x) 1 = h(x)u(x)g%[gu(xyh(x)i] = h(x)[v^(x) + d^]h{x)^ . (3.26) 

Writing 

v,(x) = iv;(x)a a = i ( <^ x) ) , v±{x) = vl(x) T (3.27) 
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and using eq.(3.17) this implies 

vl{x)' = vfc) - d^a{x), v±{x)' = exp(±2ta(x))v±(x). (3.28) 

Hence, v^(x) transforms like an Abelian gauge field for U(l) s , while v^(x) represent 
vector fields "charged" under U(l) s . For later convenience we also introduce the 
Hermitean charged vector field 

V,(x) = vfr)*! + vl(x)a 2 = v;(x)a + + v~{x)a_ = ( ^ <^ ) , (3.29) 

where a± = |(<Ti±i<7 2 ) are raising and lowering operators of spin. Under the SU(2) S 
symmetry the charged vector field transforms as 

V„{x)' = h(x)V^(x)h(xy. (3.30) 

The magnon action can also be written as 



% M ] = J d 2 x dt 2p s (v+v. + = / d 2 x dt p s Tr [u/U, + ±V t *V t 



. (3.31) 



It should be pointed out that the fields v^(x) do not represent independent degrees 
of freedom, but are composed of magnon fields. In particular, what looks like a 
mass term for a charged vector field is indeed just the kinetic term of a massless 
Nambu-Goldstone boson. 



3.4 Discrete Symmetries of Composite Fields 

Under the displacement symmetry D the composite vector field transforms as 

D v„{x) = t{x)[v^{x) + d^r(xy D ^(x) = -vfa) + d^x), 
D v±(x) = - eM^M^)K(x), D V,(x) = r(x)V,(x)r(xy. (3.32) 

Similarly, under the symmetry D' one finds 

D 'v^x) = v^x)* D 'vl{x) = -vl(x), 

D 'v±(x) = -„*(*), D 'V,(x) = -V,(x)*. (3.33) 

This is exactly how an ordinary non-Abelian gauge field behaves under charge con- 
jugation in particle physics. 

Under the 90 degrees spatial rotation O the composite field v^(x) transforms as 

°v i {x) = EijVjiOx), °v t (x) = v t (Ox), (3.34) 

while under the reflection R one obtains 

R Vl (x) = v x {Rx), R v 2 (x) = -v 2 (Rx), R v t (x) = v t (Rx). (3.35) 
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Finally, under the time-reversal symmetry T the field t> M transforms as 

T Vi(x) = D Vi(Tx), T v t {x) = - D v t {Tx) =>• 

T vf{x) = -vf(Tx) + diip(Tx), T vi{x) = vfiTx) - d t ip(Tx), 

T vf{x) = — exp(^f2i(p(Tx))v^(Tx), T vf(x) = exp(^f2i(p(Tx))vf(Tx), 

T Vi(x) = T(Tx)Vi(Tx)r(Tx)\ T V t (x) = -T{Tx)V t {Tx)T{Tx)\ (3.36) 

and under its unbroken variant T' one finds 



T 'vi{x) = 


D 'v t (Tx), 


T 'v t (x) = - D 'v t (Tx), 




T 'vf(x) = 


-vf(Tx), 


T 'vKx)=vf(Tx), 




T 'vf{x) = 


-v?(Tx), 


T 'vf{x)=vf{Tx) ) 




T 'V t (x) = 


-V t (Tx) T , 


T 'V t (x) = V t (Tx) T . 


(3.37) 



Note that the upper index T on the right denotes transpose, while on the left it 
denotes time-reversal. The above relations are equivalent to time-reversal of an 
ordinary non-Abelian gauge field. 



3.5 Alternative Representation of Magnon Fields 

We have used two equivalent representations of the magnon field in terms of the unit- 
vector e(x) and in terms of the projection matrix P(x). There is a third equivalent 
representation in terms of a complex doublet 

z{x) = (llU )' *( x y = M x ) m Mxn 

z(x)h(x) = ^(rr)! 2 + \z 2 (x)\ 2 = 1, (3.38) 

which is related to the other two representations by 

e(x) = z(x) J< az(x) =^ 
ei(x) = z 1 (x)*z 2 (x) + z 2 (x)*z 1 (x), 
e 2 (x) = i[z 2 (x)* z x (x) - z 1 (x)* z 2 (x)], 
e 3 (x) = {z^x)] 2 - \z 2 (x)\ 2 , 

u( \ ( \ ( \\ ( ki( x )| 2 z 1 (x)z 2 (x)* \ , s 

Fix) = z{x)z{xy = ' ' w I A,' . (3.39) 
v ' w w \ z 2 (x)zi(x)* \z 2 {x)\ 2 J v ' 

The field z(x) is defined in terms of e(x) or P(x) only up to a U(l) s gauge trans- 
formation 

z(x)' = exp(i(3(x))z(x). (3.40) 
It is therefore necessary to also introduce the auxiliary real-valued U(l) s gauge field 

a^x) = ^[ziptfdpzix) - d^x^ z(x)\, (3.41) 
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which under the symmetry of eq.(3.40) transforms as 

a^x)' = a ll {x)+d li p{x). (3.42) 

The complex doublet z(x) is closely related to the field u(x). Fixing the gauge 
freedom of eq.(3.40) such that z±(x) is real and non-negative, it is easy to show that 

/ zx(x) z 2 (x)* 
u{x) = \\ \ \ 

y -z 2 {x) zi(x) 

Hence, the description in terms of complex doublets z(x) and an additional auxiliary 
gauge field a M (x) is physically equivalent to what we described before. It should again 
be pointed out that a^(x) (or equivalently v^(x)) does not represent a dynamical 
Abelian gauge field, but is simply a composite field constructed from the underlying 
magnon field P(x). 



(x) 



(x). 



(3.43) 



3.6 Baby-Skyrmions 



It is interesting to note that magnon fields support topological solitons known as 
baby-Skyrmions — a lower-dimensional variant of the Skyrme soliton which repre- 
sents a baryon in the low-energy pion effective theory for QCD [83]. Baby-Skyrmions 
are solitons whose topological charge 



B — -!- J d 2 x Eije- (dicx dje), 



(3.44) 



defined at every instant in time, is an element of the homotopy group n 2 [S' 2 ] = Z. 
The corresponding topological current 



M x ) = ^ £ ^ P e{x) ■ [d u e(x) x d p e(x) 

07T 



(3.45) 



is conserved, i.e. d^j^ = 0, independent of the equations of motion. Baby-Skyrmions 
are massive excitations inaccessible to the systematic low-energy expansion of chiral 
perturbation theory. Still, the existence of the conserved current j^{x) may have 
physical consequences even for the pure magnon dynamics. 

Under the various symmetries the topological current transforms as 



SU(2) S 
SU(2) Q 
D 
D' 
O 
R 
T 

r 



D J»(x) = -j^x), 
D 'jn(x) = -jp(x), 

°3t (x)=j t (Ox), °ji (x) = Eijjj (Ox), 
R 3t(x) = -jt(Rx), R 3i(x) = -ji(Rx), R j 2 (x) = j 2 (Rx), 
T 3t(x) = -j t (Tx), T 3l (x) = 3l (Tx), T j 2 (x) = j 2 (Rx), 
T 'jt(x) = -j t (Tx), T ' 3l (x) = 3l (Tx), T 'j 2 (x) = j 2 (Tx). (3.46) 
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One might be tempted to add a term j^{x)v^{x) to the magnon Lagrangian because 
this is how an Abelian gauge field couples to a conserved current. Indeed, this term 
is invariant under SU(2) S , SU(2)q, D, D' , and O. However, it violates the reflection 
and time-reversal symmetries R, T, and T' and is hence forbidden. 

There is another non-trivial homotopy group, iyS' 2 ] = Z, which is relevant 
for baby-Skyrmions. It implies that space-time-dependent magnon fields fall into 
distinct topological classes characterized by the Hopf number H[e\ e r^fS* 2 ] = Z. 
In 2 + 1 dimensions baby-Skyrmions can be quantized as anyons characterized by 
a statistics angle 9 [84]. The cases 9 = and 9 = n correspond to bosons and 
fermions, respectively. Including the Hopf term, the magnon path integral takes the 
form 

Z = jve exp(-S[e\) exp(i9H[e\). (3.47) 

The Hopf term also changes sign under R, T, and T'. Hence, exp(iOH[e\) is invariant 
only if 9 is or n. Consequently, in an antiferromagnet with exact R, T, or T' 
symmetries baby-Skyrmions can only be quantized as bosons or fermions. For the 
antiferromagnetic quantum Heisenberg model it has been argued that no Hopf term 
is generated [85-89]. Hence, in that case the baby-Skyrmions should be bosons. 



4 The Hubbard Model in a Magnon Background 
Field 

The half-filled ground state of the Hubbard model plays a similar role as the Dirac 
sea in a relativistic quantum field theory. In particular, any fermion added to a 
half-filled state will be denoted as an electron, while any fermion removed from 
such a state represents a hole. In this section we couple a background magnon 
field to the microscopic degrees of freedom of the Hubbard model. In this way 
composite operators are constructed which transform exactly like the fields of the 
effective theory. Hence, the effective fields inherit their transformation properties 
under symmetry operations from the Hubbard model degrees of freedom. 

4.1 Fermion Operators in a Magnon Background Field 

In order to analyze the transformation properties of the electron and hole fields, as 
an intermediate step between the microscopic and effective descriptions, we first add 
a continuum magnon background field P(x) to the Hubbard model by hand. The 
corresponding diagonalizing unitary matrix field u(x) is used to turn the matrix- 
valued Hubbard model operator C x of eq.(2.24) into new operators ^ and ^ 
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defined on the even and odd sublattices, respectively 

In order to achieve a consistent representation of the underlying antiferromagnetic 
structure, it is unavoidable to explicitly split the degrees of freedom according to 
their location on sublattice A or B. In this context it may be interesting to consider 
the electron-hole representation of the Hubbard model operators discussed in ap- 
pendix A. The operators ip x ^ obey standard anticommutation relations. It should 
be noted that here the continuum field u(x) is evaluated only at discrete lattice 
points x. 

The new lattice operators inherit their transformation properties from the oper- 
ators of the Hubbard model. According to eqs.(3.16) and (2.8), under the SU(2) S 
symmetry one obtains 

= u{x)'C' x = h(x)u(x)g^gC x = h(x)^ B . (4.2) 

In components this relation takes the form 

^x±' = exp(±ia(x))tp x ^ . (4.3) 

The components ip x ± do not simply correspond to spin up and spin down with 
respect to an arbitrarily chosen global quantization axis. Instead they correspond 
to spin parallel (+) or antiparallel (— ) to the local staggered magnetization. This 
follows from considering global symmetry transformations g G U(l) s in the unbro- 
ken subgroup of SU(2) S which describe rotations around the spontaneously selected 
direction of the staggered magnetization vector. In that case, according to eq.(3.16), 
h(x) = g becomes a global transformation as well and eq.(4.3) shows that ipx± m_ 
deed has spin parallel or antiparallel to the direction of the staggered magnetization. 

Similarly, under the SU(2)q symmetry one obtains 

QqAB = Q u ( x fc x = u(x)C x n T = ^,B n T (44) 

In particular, under the U(1)q subgroup of SU(2)q the components transform as 

Q ^x± = expM^f- (4.5) 
Under the displacement symmetry the new operators transform as 

D ^ B = D u(x + i)C x+l a 3 = r(x + i)u{x + i)C x+l a 3 = r(x + *)*f> 3 , (4.6) 



X 
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where t(x) is the field introduced in eq.(3.21). Expressed in components this implies 

= TeMTMx + <))^ >T . (4-7) 

Similarly, under the symmetry D' one finds 

D '^ B = D 'u(x + i)(ia 2 )C x+l a 3 = u(x + i)*(ia 2 )C x+ .a 3 = («7 2 )tff£<7 3 . (4.8) 

Here we have used u(x + i)*(io" 2 ) = {ia2)u{x + i). Again, expressed in components 
this relation takes the form 

= ±€i T - (4-9) 

We have seen before that the symmetry D' acts on the composite field v^(x) exactly 
like charge conjugation in particle physics. However, it should be noted that D' acts 
on the electron and hole fields in a different way than the usual charge conjugation 
of a relativistic Dirac fermion which interchanges electrons and positrons. In partic- 
ular, D' does not interchange electrons and holes. Instead, it flips the spin of both 
electrons and holes from + to — and vice versa. Indeed, the spin is the "charge" 
that couples to the composite gauge field of eq.(3.25) constructed from the magnon 
field. 

In the condensed matter literature on high-temperature superconductivity the 
concept of spin-charge separation (whose existence is established for some systems 
in one spatial dimension) has often been invoked. The idea is that there may be 
quasi-particles — so-called holons — which carry charge but no spin, as well as 
so-called spinons which are neutral and carry spin \. In order to avoid confusion 
between holons and the holes of our effective theory, we like to make a few comments: 
one might think that the fermion operator ^ x does not carry spin since it does 
not transform with the global spin transformation g G SU(2) S . However, the spin 
symmetry is non-linearly realized and hence the fermion operator transforms with 
the local h(x) G U(l) s . Consequently, ^ X ' B still carries spin and hence does not 
represent a holon. It should also be pointed out that in the weakly coupled effective 
theory of magnons and holes there are no linearly confining forces that could form 
a spinless holon out of ^ X ' B and the magnon field z(x) of eq.(3.38). 



4.2 Formal Continuum Limit of the Hubbard Model in a 
Magnon Background Field 

In terms of the new operators the Hubbard model Hamiltonian takes the form 

x£B,i 
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-f E Tr[^f ^ 3 ] - f E n[*?tt?<7 3 ], (4.10) 

where we have introduced the parallel transporter 

V x ,i = u(x)u(x + i)t e SU(2) a , (4.11) 
which transforms under SU(2) S as 

^,i = M^)V,,iMa; + <)t. (4.12) 
For smooth magnon fields we can put 

u(x) = u(x +'-)- ^d iU (x +'-) + jdM* +\) + 0(a 3 ), 

u(x + i) = u(x + % -) + ^d iU (x + % -) + jd*u(x + \) + 0(a 3 ), (4.13) 

where a is the lattice spacing. Similar expressions hold for the other fields. Using 
the unitarity of u(x + |) one can show that the lattice parallel transporter reduces 
to 

V Xti = 1 + avi(x + l -) + ^(x + l -f + £>(a 3 ), (4.14) 

with Vi(x) given by eq.(3.25). Note that both the continuum field Vi(x) and the 
lattice parallel transporter field V x ^ transform locally only with the unbroken U(l) s 
subgroup and not with the full SU(2) S symmetry. 

In the continuum limit we make the replacements 

£> E - A/^, < B - V2a* A > B (x). (4.15) 

x(zA x&B 

The factor | in front of the integral accounts for the fact that each sublattice covers 
only half of the space. Similarly the factor \/2a in the definition of the continuum 
field ty A ' B (x) arises because there is only one degree of freedom of a given type 
A or B per area 2a 2 . The components ip± B (x) of ^/ A ' B (x) again obey standard 
anticommutation relations, however, with the Dirac ^-function of the continuum 
theory instead of the Kronecker (^-function of the lattice. It should be noted that, 
due to the antiferromagnetic order, the number of degrees of freedom per continuum 
point is twice as large as the number per lattice point. Taking the formal continuum 
limit a — > (and ignoring an irrelevant constant) the Hamiltonian of eq.(4.10) takes 
the form 

H = J d 2 x {MTr[^ At ^ B ] + ^Tr[A^ At A^ B ] 

+iKTr[D i ^V i ^ B + D^V^] + 7VTr[* A V^ B ] 
+ 9jj t \^a\^a^a^a + ^b^b^b^b^ _ ^Tr[^ At ^"V 3 + ^ B ^ B a 3 ]}. 

(4.16) 
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It should be noted that, due to the structure of ^ A ' B (x), the individual terms are 
Hermitean. In the above expression V^x) is the field defined in eq.(3.29) and the 
covariant derivatives are given by 

D^ B {x) = (d, + ivl(x)a 3 )V A > B (x), 

D^ A ' B \x) = [D^ A > B {x)\< = d^ A ' B \x) - ^ A ' B \x)ivl(x)a z . (4.17) 

In terms of the fundamental parameters t and U and the lattice spacing a of the 
Hubbard model one obtains 

M = -4t, M' = ^—, K = ta 2 , N = ta 2 , G = 2Ua 2 . (4.18) 
2ta 2 

It should be noted that (in contrast to a relativistic theory) the kinetic mass M' 
is in general different from the rest mass M. The Hamiltonian from above resem- 
bles some (but not all) terms in the action of the effective theory to be constructed 
below. However, the coupling constants resulting from the formal continuum limit 
get renormalized and will hence be replaced by a priori unknown low-energy pa- 
rameters in the effective action. The values of the low-energy parameters can be 
determined in experiments with cuprate materials or through numerical simulations 
of a microscopic Hubbard-type model. 



5 Effective Theory for Magnons and Charge Car- 
riers 

The low-energy effective theory for magnons is analogous to chiral perturbation the- 
ory for pions in QCD. In QCD the baryon number B is a conserved quantity. Thus 
one can investigate the low-energy QCD dynamics separately in each baryon num- 
ber sector. Ordinary chiral perturbation theory operates in the B = sector. The 
low-energy physics in the B = 1 sector involves a single nucleon interacting with soft 
pions. The low-energy effective theory describing these dynamics is known as baryon 
chiral perturbation theory [69-73] . Similar effective theories have been constructed 
for the B = 2 [52, 53] and B = 3 sectors [55, 60] in the context of nuclear physics. 
Even nuclear matter (i.e. a system with non-zero baryon density) has been studied 
with effective theories [62-66]. The condensed matter analog of baryon number is 
electron (or hole) number (or equivalently electric charge) which is obviously also 
conserved. In analogy to QCD it is hence possible to construct a low-energy effec- 
tive theory describing the interactions of soft magnons with charge carriers doped 
into an antiferromagnet. Most high-T c materials result by hole-doping of quantum 
antiferromagnets, but the effective theory also applies to electron-doping. The key 
observation is that the spontaneously broken SU(2) S spin symmetry is non-linearly 
realized on the electron or hole fields and appears as a local U(l) s symmetry in the 
unbroken subgroup. This is analogous to baryon chiral perturbation theory in which 
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the spontaneously broken SU (2)l®SU (2)r chiral symmetry of QCD is implemented 
on the nucleon fields as a local SU{2) L=R transformation in the unbroken isospin 
subgroup. 



5.1 Effective Fields for Charge Carriers 



In the low-energy effective theory we will use a Euclidean path integral description 
instead of the Hamiltonian description used in the Hubbard model. Consequently, 
the Hermitean conjugate lattice operators ip x ± are then replaced by Grassmann 
numbers ip±' B \x) which are completely independent of ip±' B (x). Therefore, in the 
effective theory the electron and hole fields are represented by eight independent 
Grassmann numbers ip± B (x) and ^± B \x) which can be combined to 



ip A (x) ip_ (x) 




—ip B {x) 



X 



(5.1) 



In order to avoid confusion with relativistic theories, we do not denote the conjugate 
(x). For notational convenience we also introduce the fields 



fields by ip ± 
V A \x) 



i>-(x) -4> A {x) 



m B \x) 



^+\ x ) 



ip B (x) 



-ip_(x) ^ B {x) 



(5.2) 



It should be noted that ^/ A ' B ^(x) is not an independent field, but consists of the 
same Grassmann fields ip± B (x) and ip± B \x) as ty A ' B (x). 

It should be pointed out that, since they emerge dynamically, the continuum 
fields of the low-energy effective theory can not be derived explicitly from the lattice 
operators of the microscopic Hubbard model. Still, the Grassmann fields ty A ' B (x) 
describing electrons and holes in the low-energy effective theory transform just like 
the lattice operators ^ A,B discussed before. In contrast to the lattice operators, 
the fields ^/ A ' B (x) are defined in the continuum. Hence, under the displacement 
symmetries D and D' one no longer distinguishes between the points x and x + i. 
As a result, the transformation rules of the various symmetries take the form 



SU(2) S 
SU(2) Q 
D 
D' 



^ B (x)' - 
^ A > B (x) 
D ^ B (x) 
D '* A ' B (x) 



h(x)V A > B (x), 

= ^ B (x)n T , 

- t(x)* b ' a (x)ct 3 , 
= (ia 2 )V B > A (x)a 3 . 



(5.3) 



In components the symmetry transformations read 



SU(2) S 
U(1)q 
D 
D' 



ip A ' B (x)' = exp(±ia(x))ip A ' B (x) , 
Q ip± B {x) = exp(iu)ip± B (x) , 
D ipi' B (x) = =fexp(^i(p(x))^' A ( x )> 
D '^ B (x) = ±^\x). 



(5.4) 
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Under the space-time symmetries, i.e. under the 90 degrees rotation O, the reflection 
R, time-reversal T, and its unbroken variant T' the fermion fields transform as 



: 


°V A ' B (x) = 


V A ' B (Ox), 


R : 


R q? A > B (x) = 


V A > B (Rx), 


T : 


T ^ B {x) = 


T(Tx)(ia 2 )* A ^(Tx) T a 3 , 




T m A ' B \ x ) = 


-- -a 3 ^ A ' B (Tx) T (ta 2 )^r(Txy 


T' : 


T '^ B (x) = 


-V A ^(Tx) T a 3 , 




T '^ B \x) - 


= a 3 ^ B (Tx) T . 



Again an upper index T on the right denotes transpose, while on the left it denotes 
time-reversal. The form of the time-reversal symmetry T in the effective theory 
with non-linearly realized SU(2) S symmetry follows from the usual form of time- 
reversal in the Euclidean path integral of a non-relativistic theory in which the spin 
symmetry is linearly realized. The fermion fields in the two formulations just differ 
by a factor u(x). In components the previous relations take the form 



: 


°^ B {x) = 


^ B (Ox), 




R : 


R ^ B (x) = 


^ B (Rx), 




T : 


^ A ' B (x) = 
T i>t' B \x) = 


exp(^(Tx))^' B \Tx), 
-- -exp(±^(Tx))^±' B (Tx), 




T' : 


T '^' B (x) = 


-^> B \Tx), 






T '^ B \x) - 


= ^ B (Tx). 


(5.6) 



It should be noted that the components + and — (denoting spin parallel and an- 
tiparallel to the direction of the staggered magnetization) are not interchanged under 
time-reversal. While both the spin of the fermion and the staggered magnetization 
change sign under time- reversal, the projection of one onto the other does not. 

The action to be constructed in the next section must be invariant under the 
internal symmetries SU(2) S , U(1)q (or even SU(2)q), D and D' , as well as un- 
der space-time translations and the other space-time symmetries O, R, and T (or 
equivalently T'). 

The fundamental forces underlying condensed matter physics are Poincare- inva- 
riant. However, some of the space-time symmetries may be spontaneously broken 
by the formation of a crystal lattice. The resulting Nambu-Goldstone bosons are 
the phonons, which play a central role in ordinary low-temperature superconductors 
by providing the force that binds Cooper pairs. In high-T c superconductors, on the 
other hand, it is expected that phonons alone cannot provide the mechanism for 
Cooper pair formation. In the Hubbard model (and also in our effective theory) 
phonons are explicitly excluded because one imposes a rigid lattice by hand. This 
does not only break continuous translations and rotations down to their discrete 
counterparts; it also breaks space-time rotations. In a relativistic context these 
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would be the boosts of the Poincare group. In a non-relativistic theory the lattice 
explicitly breaks Galilean boost invariance, thus providing a preferred rest frame (a 
condensed matter "ether"). As a consequence, the magnon- mediated forces between 
a pair of electrons or holes may depend on the center of mass momentum of the pair. 
In the actual high-T c materials Galilean (or more precisely Poincare symmetry) is 
spontaneously (and not explicitly) broken. If phonons play an important role in 
the understanding of high-temperature superconductivity, one should construct an 
effective theory of spontaneously broken (and thus non-linearly realized) SU (2) s and 
Galilean symmetry which would automatically include both magnons and phonons. 
This is indeed possible and presently under investigation using the techniques of low- 
energy effective field theory. In the present paper we assume that phonons play no 
major role in the cuprates. In that case, it is legitimate to break Galilean invariance 
explicitly instead of spontaneously. 



5.2 Effective Action for Magnons and Charge Carriers 

We now construct the leading terms in the effective action of magnons and electrons 
or holes. The effective theory provides a systematic low-energy expansion organized 
according to the number of derivatives in the terms of the effective action. We 
decompose the effective Lagrangian into an S77(2)Q-invariant part £ and an SU(2)q- 
breaking (but still £7(1 ^-invariant) part C The contributions C nunun ^ and C nununi> 
to the effective Lagrangian are classified according to the number of time- derivatives 
n t , the number of spatial derivatives n i; and the number of fermion fields they 
contain. The total action is then given by 

S[lp± B \lp± B ,P}= d 2 X dt {£>n t , ni ,n^ +£n t ,n i ,nj (5.7) 

and the partition function takes the form 

Z = J V^V^i' B VPexp(-S[^ B \^± B ,P])- (5.8) 

Until now we have constructed the effective action in the Q = sector, i.e. for 
a half-filled system which is described entirely in terms of magnons. Since antifer- 
romagnetic magnons have a "relativistic" dispersion relation (with the spin-wave 
velocity c playing the role of the velocity of light), in pure magnon chiral perturba- 
tion theory one counts temporal and spatial derivatives as being of the same order. 
The leading contributions of eq.(3.5) take the form 

£2,0,0 = ^Tr[d t Pd t P], £0,2,0 = PsTrfoPdiP]. (5.9) 

Next we consider terms quadratic in the fermion fields. These contribute to the 
scattering of magnons off electrons or holes in the Q — ±1 sectors and they generally 
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describe the propagation of charge carriers in an antiferromagnet with \Q\ > 1. 
In contrast to magnons, electrons or holes are massive and have a non-relativistic 
dispersion relation. Hence, it is natural to count one temporal and two spatial 
derivatives as being of the same order. In order to count derivatives consistently, in 
the Q 7^ sectors it may thus be necessary to also consider the pure magnon term 
£2,0,0 with two temporal derivatives as being of higher order. The leading order 
terms without any derivatives which are Hermitean and invariant under SU(2) S , 
SU(2)q, D, and D' as well as under the space-time symmetries O, R, T, and T' 
take the form 

£ 0)0j2 = MiTr[tf At tf B ] + ^Tr^Vs^ 4 - tf st <r 3 tf B ] 

= m x V+ + i>-Hl + i>lHi + i>- H- ) 

+M 2 (^Vi - V'-V- - + i>- V-)- (5-io) 

The mass parameters M 1 and M 2 (as well as all other low-energy parameters to be 
introduced below) take real values in order to ensure Hermiticity of the correspond- 
ing Hamiltonian. It should be noted that 

Tt[^ a ^ a ] = Tt[^ b ^ b ] = 0, 

Tr[* At * B ] = Tr[V B ** A ], (5.11) 

due to the anticommutativity of the Grassmann fields. When we impose only the 
generic U(1)q but not the full SU(2)q symmetry, one more fermion mass term can 
be added 

= m(Vf^ + V-V_ + (5.12) 

This term can be absorbed into a redefinition of the chemical potential. Remarkably, 
no other fermion mass terms (consistent with the SU(2) S , U(1)q, D, D', T, and T' 
symmetries) exist. In particular, it is useful to note that 

Tr[# At <7 3 * A <7 3 ] = Tr[^V 3 ^ B CT 3 ] = 0. (5.13) 

The terms with one temporal derivative are given by 
£1,0,2 = ±Tr[* A tD t * A + * B iD t * B ] 

2 

= itfDtrft + V- f AV>- + j?D t rft + ^ D t ^ B _ 
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Here V t is the field defined in eq.(3.29) and the covariant derivatives are those of 
eq.(4.17). In components they take the form 

D li ^' B (x) = (d li ±ivl(x))ipi' B (x), 

D^ B \x) = (d, T ivl(x))i>£> B \x). (5.15) 

Note that vf as well as vf (and hence V t ) count like one temporal derivative because 
these composite fields indeed contain one time-derivative of the magnon field. 

When one derives the Euclidean path integral from the Hamiltonian formulation 
of the effective theory, the term V'+^9tV'++V'-^tV'-+V'+^tV'++V'-^tV'- arises from 
the pairs of anticommuting fermion operators. It should be noted that there are two 
more S'?7(2)Q-breaking but £/(l)g-invariant terms with a single time- derivative 

^Tr[^a 3 D t ^ A a 3 - ^a 3 D^ B a 3 \ 

= itfDttf - ^ f AV>- - ^D^l + ^ D t ^ B _ , 
±Tr[* A 1 D t * B a s + ^ D t ^ A a 3 ] 

= ^D t 4> B + ip A ^D t ip B + ip B ^ Dfip A + ^D t 4> A . (5.16) 

These terms need not be included in the effective Lagrangian, since they would not 
imply canonical anticommutation relations in the Hamiltonian formulation. In any 
discussed in appendix B, if one does include these terms they can again be 
removed by an appropriate field redefinition. 

Interestingly, there is only one more term that violates the SU(2)q symmetry 
but still respects the U(1)q symmetry 

A,o,2 = XTr[^V t ^ B a 3 ] 

= X(^v^ B + ^v^ + ^v^ A + ^v^ B ). (5.17) 

Further potential contributions are absent because, for example, 

Tr[* At A^V 3 ] = Tr[^ m D t ^ B a 3 ] = 0, 
Tr[^ A ^V t ^ A a 3 ] = Tr[^ B ^V t ^ B a 3 ] = 0, 

Tr[* At <T 3 V t V A <T 3 ] = Tr[^a 3 V t ^ B a 3 ] = 0. (5.18) 



Terms with a single spatial derivative are forbidden due to the reflection symme- 
try R and the 90 degrees rotation symmetry O of the quadratic spatial lattice of the 
underlying microscopic system. The terms with two spatial derivatives are given by 

£0,2,2 = ^Tr[Di^ D^ B ] + -^^[D^a 3 D^ A - D^a 3 D^ B ] 
+iK 1 Tr[D i y A *V i y B + D^V^] 
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+iK 2 Tr[D i V A i<T 3 W A - DiV^azW 3 ] 

+N 1 Trl^ViV i y B ] + — Ti^Vi^V^ - ^ Bf V> 3 V^ B ] 

2 

= ^{D^D^l + D^D^i + D^D^ B + D^D^ A ) 
+ ^(A^f A^ - D^D^ A - D^D^l + D^D4 B ) 

+iK 1 {Dnl>?vf<l>- ~ V>- VA^ + AV>-VV>+ - VA^- 
+D l ^vf^ A - ^v-D^l + A^-V #f - #JV Wf) 

+tK 2 (D t ^v^ A - ^v~D^ A - AV>- + ^v+D^ A 
-D^vf^ B + ^v'D^l + D^v'^l - ^vfD^ B ) 

+N 1 (^v+v~^ B + ^v~vt^ B + ^vfv~^ A + ^- V^_) 

-N 2 (^vtv~^ A - ^v~vf^ A - ^vfv-^l + V- V^-)<5.19) 

Note that the imaginary unit i in front of the terms proportional to K\ and is 
necessary to ensure that the corresponding Hamiltonian is Hermitean. In principle, 
terms containing AA and A^i could also be written down. However, upon partial 
integration, up to irrelevant surface terms they lead to the same Euclidean action as 
the terms constructed here. Since the doped electrons or holes are non-relativistic, 
there is no reason why the kinetic mass parameters M[ and M' 2 should agree with 
the rest mass parameters M\ and M 2 . In addition, there are again terms that break 
the SU{2)q symmetry but leave the U(1)q symmetry intact 

£0,2,2 = -^TrfA^ 41 A^Vs + A^A^s] 
Am' 

+tn 1 Tr[D l ^a 3 V^ B a 3 + A* Bt K<7 3 * A <7 3 ] 
+tK 2 Tr[D t ^V t V A cr 3 + D^V^ B a 3 ] 

+^Tr[^V^ A a 3 + m^ViV^az] 

= ^7 (D t ^D^ A + D^D^ A + D^l f + A^_ 1 ' A V'- ) 

+i«i(A^V^_ - V^VA^ - D^v'^l + ^VA^- 
- AV>+ + ^v-D^l + D^v-^ A - ^vfD^ B ) 

+iK 2 {D^vti> A - ^v-D^i + A^- V#J - ^VA^- 
+A^+ - ^v'D^l + A^VV>+ - ^v+D^ B ) 

+v(^v+v-^ A + V- + 1>?v+vrifi + ^ B Jv-vti> B ). (5.20) 

Next we consider terms quartic in the fermion fields which describe short-range 
interactions between the charge carriers. To lowest order there are five linearly 
independent 4-fermion contact interaction terms 

C 4 = ^lTr[^ A ^ A ^ A ^ A + y B ^ B ^ B ^ B ] + ^Tr[* At * A * Bt * B ] 



31 



+G 5 Tr[^ A V 3 ^ A ^ At ^ B - ^a 3 ^ B ^f B ^ A ] 

(5.21) 

It is interesting to note that 

Tr[$! A ^ B ty B ^ A + iIr A t^ijfBtijfB _|_ 2^ A tijr B ^ A t^B] = q ? 
Tr[^ A t^^ A + ^t^^t^A + 2 ^ A ^a 3 ^ A ^ B ^a 3 ^ B ] = 0, 
Tr[^ At ^ B ^ Bt ^ A - ^ a ^ a ^j b ^ b ] = 2(Tr[^ At ^ B ]) 2 . (5.22) 

Together with other relations similar to these ones, this implies that the terms listed 
above form a maximal linearly independent set. 

Again, there are additional terms that are invariant under U(1)q but not under 
SU(2) Q 

£ 0A4 = ^Tr[^ At a 3 V A V B *V B a 3 - ^ m a 3 ^ B ^ A ^ A a 3 ] + ^Tr[^ A ^ A a 3 ^ B ^ B a 3 ] 
+# 3 Tr[^ At ^ A o-3^ At ^ + ^ B ^ B a 3 ^ B ^ A } 

(5.23) 

One may note that, for example, 

Tr[^ At ^ A ^ Bt ^ B + ^ A ^ A a 3 ^ B ^ B a 3 + 2^ A ^ B a 3 ^ B ^ A a 3 ] = 0. (5.24) 

Together with further relations of a similar kind, this implies that there are no other 
linearly independent 4-fermion terms that obey the relevant symmetries. 



32 



For completeness, let us also construct the terms containing six fermion fields 
and no derivatives. The SU (2 ^-invariant 6-fermion terms can be written as 

£0,0,6 = ^Tr[* A t<73^ At <73^ B V 3 * S - ^V 3 ^ B V 3 ^ A V 3 ^ A ] 
+ ^Tr[^ j4t ^ B ^ At ^ B ^ At ^ B ] 

= # 1 (^V^+V+V'-V- + V B 

V-V_^iV+ + ^V^+V+^-V.)- (5.25) 

It is interesting to note that 

Tr[* A ia 3 * A * A ia 3 * A * B *<T 3 * B - ^a 3 ^ Bjf a 3 ^ B ^a 3 ^ A ] 

= -Tr[^a 3 ^ A - ^V 3 ^ B ]Tr[^ At a 3 ^ B V 3 ^ B ]. (5.26) 

In addition, there is one S*?7(2)Q-breaking (but [/(l)g-invariant) 6-fermion term 
£0,0,6 = ^Tr[^ A ^ A a 3 ^ A ^ A a 3 ^ B ^ B a 3 + * b ** b <t 3 * b ** b <t 3 * a *V a <t 3 ] 

V_ + V-)- (5.27) 

Finally, the only 8-fermion term with no derivatives takes the form 
£ 0A 8 = ^ T T [^ A ^ A ^ B ^ B ^ A ^ A ^ B ^ B ] 

= -1 viV^+V+V'-V-V'-V-, (5.28) 

which is S , f/(2)Q-invariant. It may be noted that 

Tr[* At * A * Bt * B * At * A * st * B ] + i(Tr[^ Af ^ A \l/ Bt \l/ B ]) 2 = 0. (5.29) 

2 

No 5*^7 (2)Q-breaking 8-fermion term without derivatives exists, such that £0,0,8 = 0. 
Terms with more than eight fermion fields vanish due to the Pauli principle, unless 
one includes derivatives. Since such terms are of higher order than those without 
derivatives, they will not be constructed here. When one wants to address questions 
for which the short-distance forces between charge carriers are essential, it will be 
necessary to consider such terms. While constructing them is a straightforward 
exercise, it is not very illuminating and will hence be omitted at this stage. 

The 4-, 6-, and 8-fermion contact terms parameterize short distance interactions 
with a large number of undetermined low-energy constants. Since this limits the pre- 
dictive power of the effective theory at short distances, it is natural to concentrate 
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on long-distance forces between the charge carriers. For example, one-magnon ex- 
change mediates a long-range force that is unambiguously predicted by the effective 
theory in terms of just a few low-energy parameters. 

It should be mentioned that there are many equivalent ways of rewriting the 
various contributions to the action in terms of traces. Hence, the above choices of 
terms are to some extent arbitrary. It is important that the selected terms form 
a maximal linearly independent set. For example, all determinants or products of 
traces of fermion fields can be written as linear combinations of the traces listed 
above. To verify the completeness and linear independence of the selected terms is a 
non-trivial task which was addressed by extensive use of the algebraic manipulation 
program FORM [90]. 

It is straightforward to include the fermion chemical potential \i in the effective 
theory. It appears as the temporal component of a purely imaginary U(1)q gauge 
field and thus manifests itself in an additional contribution to the covariant derivative 

D t * A > B {x) = d t ^ A > B (x) + ivl(x)a^ A > B (x) - iiV a > b (x)<j 3 . (5.30) 

Before one can do consistent loop-calculations in the low-energy effective theory at 
non-zero Q or at non-zero /i one must develop a power-counting scheme, e.g. along 
the lines of [73]. This will be the subject of a future publication. 



5.3 Dispersion Relations of Electrons and Holes 

In this subsection, as an application of the effective theory, we consider the dispersion 
relations of the charge carriers. For this purpose we switch off the magnon field (i.e. 
P(x) = |(1 + <7 3 ) =>- u(x) = l,f At (x) = 0) and consider the propagation of free 
charge carriers in the antiferromagnetic medium. In the absence of SU (2)Q-breaking 
terms, the Lagrangian (quadratic in the fermion fields) then reduces to 

C = ^d t ^ A + ^d t ^ A + ^d^ B + ^d^ B 




Mi 




-M 2 


H 


i 

2M' 2 
1 


1 

2M[ 
1 


2M[ 


2M^ 


1 

2M' 
1 

2M[ 


1 

2M[ 
2W 2 




+(^U^t) 2 ^ ^ °Tb • (5-31) 



The eigenstates of free particles propagating with a 2-d momentum vector p arise 
as the eigenvectors of the matrices 

H+{P) = 1 M 1 + 4 -M 2 -A 

1V±1 ^ 2M[ 1V12 2ML, 
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Due to the lack of Galilean invariance the eigenvectors depend on p 2 , i.e. the prob- 
ability for an electron or hole to be found on the A or B sublattice depends on 
the momentum. As a consequence of the displacement symmetries D and D' the 
eigenvalues of H + {p 2 ) and H_{p 2 ) are the same. Both matrices have two eigenvalues 




(5.33) 



The positive energy states correspond to electrons, while the negative energy states 
correspond to holes. Not surprisingly, due to the SU(2)q symmetry electrons and 
holes have the same dispersion relation. The rest mass M and the kinetic mass M' 
are given by 

M = v /Mf + Mf, ^ = § + ^. (5.34) 
v 1 1 M' M[ M' 2 v ' 

Next we take into account the additional terms that reduce the SU(2)q symmetry 
to the U(1)q symmetry. Then there are additional contributions to the energy 

H + {p 2 )=H^p 2 )=( m + ^ ° p2 ). (5.35) 

and the corresponding eigenvalues now take the form 

Still, the energies in the + and — sectors are the same. However, the electron and 
hole dispersion relations now differ. 

At this point, we have constructed eigenstates of the free Hamiltonian with defi- 
nite continuum momentum and with definite spin projection on the direction of the 
staggered magnetization. However, unlike the eigenstates of the underlying micro- 
scopic Hamiltonian, the states of the effective theory do not have a definite lattice 
momentum. Still, the low-energy effective theory defined in the continuum knows 
about the underlying lattice structure through the realization of the displacement 
symmetries D and D' . Since the symmetry D is spontaneously broken, neither the 
vacuum nor the single particle states are eigenstates of D. Operating twice with D 
acts trivially on the fields, i.e. DD P(x) = P(x), DD ^/j_' B = and hence does 

not reveal any useful information. It is more useful to operate with the unbroken 
displacement symmetry D' . In particular, the vacuum state P(x) = |(1 + 03) is 
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invariant under D' . Still, in the way we constructed them, the electron or hole 
states of the effective theory are not eigenstates of D'. However, since states with 
spin parallel and antiparallel to the staggered magnetization are degenerate with 
each other, one can form appropriate linear combinations that are eigenstates of the 
displacement symmetry D' . Applying D' twice one obtains 

D ' D '^ B (x) = ± D '^ A (x) = -^ B (x), (5.37) 

which implies that the corresponding eigenvalue A = exp(ika) of D' obeys 

A 2 = exp(2ika) = -1 =>• ka = (5.38) 

This is reminiscent of the result, mentioned in the introduction, that low-energy 
hole states are located at lattice momenta (±|,±|) [6,7,17,34,37-40]. However, 
the comparison with these findings is subtle. In particular, the results of the exact 
diagonalization study on small [17] and of the Monte Carlo study on larger volumes 
[34] must be interpreted carefully. In a finite volume (with periodic boundary condi- 
tions), in analogy to QCD [91], both the SU (2) s spin symmetry and the displacement 
symmetry D are restored and the staggered magnetization acts as a quantum rotor 
[48]. As a result, in contrast to the infinite volume limit, the single particle states 
in a finite volume can be constructed as eigenstates of D. It is interesting to note 
that the finite volume effects that lead to the restoration of the spontaneously bro- 
ken symmetries SU(2) S and D can be understood in the framework of the effective 
theory. This requires a nonperturbative quantum mechanical treatment along the 
lines of [48,91]. 



6 Systems with Holes only 

In this section we consider the t-J model as well as its low-energy effective the- 
ory. In the t-J model holes are the only charge carriers which leads to substantial 
simplifications in the effective theory. 

6.1 The t-J Model 

The t-J model is defined by the Hamilton operator 

h = p { -* E(4<wi + i +fx ) + j£4 • s x+ . - ^5>* - !) 1 p > (e- 1 ) 

^ x,i x,i x J 

with 

c x — { ^ ] , S x — c x —c x , n x = c x c x . (6-2) 
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In contrast to the Hubbard model, in the t-J model the operators act in a restricted 
Hilbert space of empty or at most singly occupied sites. In particular, states with 
doubly occupied sites are exiled from the physical Hilbert space by the projection 
operator P. Hence, by definition, the t-J model does not allow the addition of 
electrons to a half-filled state. Consequently, the only charge carriers are holes. 

It is straightforward to show that the t-J model has the same symmetries as 
the Hubbard model. The only exception is the SU(2)q symmetry which relates 
electrons to holes in the Hubbard model, and which is absent in the t-J model. 
Still, the Abelian fermion number symmetry U{1)q remains exact in the t-J model. 



6.2 Effective Theory for Magnons and Holes 

Since, up to the SU(2)q symmetry, the t-J model has the same symmetries as 
the Hubbard model, the effective theory of the previous section also applies in this 
case. Of course, the values of the low-energy parameters will be different than 
for the Hubbard model. Still, the absence of electrons beyond half-filling leads to 
drastic simplifications. In particular, in the effective theory the absence of electrons 
manifests itself by an infinite electron rest mass. Consequently, with a finite amount 
of energy these excitations cannot be generated. As discussed in the previous section, 
the diagonalization of the mass matrices of electrons and holes yields 



m±M 2 Mi \ t _ / m ± y/ M? + M$ 

Mi mTM 2 ) ± ~y o mT sj Ml + Mi 




U ± = ( v Y , X,YeR. (6.3) 



The eigenvectors corresponding to the eigenvalue m + JMl + Mf describe electrons, 



while the ones corresponding to m — J Mf + Mf describe holes. When the electron 



rest mass m + y Mf + Mf goes to infinity, the corresponding eigenvector fields 

+ Y^l(x) = 0, Y^(x) + Xip*(x) = 0, (6.4) 

which describe electrons, must be put to zero. The orthogonal combinations 

$ + {x) = -Y^(x) + X^(x), V-(z) = X^(x) - Y^ B (x), (6.5) 

describe holes and must be kept. As a result, the number of degrees of freedom 
is reduced by a factor of two. In complete analogy to the discussion in appendix 
B one can show that the hole field ip±(x) transforms as follows under the various 
symmetry operations 

SU(2) a : V±W = exp(±ta(x))ip ± (x), 
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D 
D' 
O 
R 
T 

r 



Q 



D' 



exp(iu;)t/j±(x)i 

tp±(x) = ±ip T (x), 
°il>±{x)=MOx), 

T ij ± (x) = exp(Ti<p(Tx))iJ)l(Tx), 
T ipi{x) = -exp(±i(p(Tx))if>±(Tx), 
T 'i, ± {x) = -A{Tx), 
T '^ ± (x)=^ ± (Tx). 



(6.6) 



Hence, except for the SU(2)q symmetry, all symmetries can also be implemented 
on the hole fields alone. It should be noted that the transformation laws for tp±(x) 
result from those for ?p± B (x) simply by dropping the sublattice indices A and B. 

The absence of electron fields also drastically reduces the number of terms one 
can write down in the low-energy effective theory. In particular, the leading terms 
in the effective action now take the form 

S[i(>t, V>±, P] = J d 2 x dt {psTridiPdiP + ^d t Pd t P] + M(^J>+ + 

+N(^ + vtv~^ + + ^_v~vt^-) + G^iV+V'-^-}- (6-7) 

This form of the effective action is similar to (but not identical with) the ones of 
[11, 13,25,32,36]. In particular, in some of those works spin-charge separation was 
invoked and spinless fermions were considered. Also the role of the sublattice indices 
(which have at this stage disappeared from our description) is different in those 
approaches. Furthermore, the dynamical role attributed to the composite gauge 
field in some of those works is different than in our effective theory. It should be 
pointed out that the above effective Lagrangian correctly describes the low-energy 
dynamics of holes only if electrons are completely absent beyond half-filling (as 
it is indeed the case in the t-J model). Otherwise the general effective theory of 
the previous section with a larger number of low-energy constants (and thus with 
somewhat reduced predictive power) must be employed. 



7 Coupling to External Electromagnetic Fields 

In the following sections we will couple both microscopic and effective theories for 
antiferromagnets to external electromagnetic fields. For this purpose, we will make 
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use of an observation by Frohlich and Studer concerning the Pauli equation [78] . 



7.1 Local SU(2) S Symmetry of the Pauli Equation 

Up to corrections of order 1/M| (where M e is the electron mass) the Pauli equation 
(i.e. the non-relativistic reduction of the Dirac equation to its upper components) 
takes the form 

H9, - » + ■ E + iJL.3 ■ <?)* = --Ltf + iel- i^-S x S)>». (7.1) 

Here ty(x) is a 2-component Pauli spinor at the space-time point x — (x,t), a are the 
Pauli matrices, <&(#) and A(x) are the electromagnetic scalar and vector potentials, 
and 

E(x) = -V$(x) - d t A(x), B(x) = V x A(x), (7.2) 

are the usual electromagnetic field strengths. The first two terms on the left-hand 
side of eq.(7.1) form the U(1)q covariant derivative familiar from QED. The third 
(Darwin) and fourth (Zeeman) term on the left-hand side represent relativistic cor- 
rections. The first two terms on the right-hand side again form an ordinary U(1)q 
covariant derivative, while the third term represents the relativistic spin-orbit cou- 
pling. The Pauli equation transforms covariantly under U(1)q gauge transforma- 
tions 

G*(a;) = exp(iu(x))V(x), Q $(x) = $(x) + -d t u(x), Q A(x) = A(x) - -Wu(x). 

e e 

(7.3) 

Obviously, it is also covariant under global spatial rotations 

°^(x,t) = gV(Ox,t), °$(x,t) = $(Ox,t), °A(x,t) = T A(Ox,t). (7.4) 
Here O is a general orthogonal 3x3 rotation matrix with 

O t 3 = g^ffg, (7.5) 
where g G SU(2) S represents the rotation O G SO (3) in spinor space. 

Frohlich and Studer noticed that the Pauli equation has a hidden local SU(2) S 
spin symmetry. This symmetry becomes manifest when one writes 

iD t m = -AA*, (7.6) 

2M e ' v ; 

with the SU(2) S <g> U(1)q covariant derivative given by 

D^d^ + W^+ieA^x). (7.7) 
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The components of the non-Abelian vector potential 

W ll (x)=xWZ(x)?±, (7.8) 

can be identified as the electromagnetic field strengths E(x) and B(x), i.e. 

W t a (x) = fi e B a (x), W t a (x) = ^e tab E\x). (7.9) 

The anomalous magnetic moment /i e = g e e/2M e of the electron (where, up to QED 
corrections, g e = 2) appears as a non-Abelian gauge coupling. The Abelian vec- 
tor potential A^(x) is the usual one, except for a small contribution to the scalar 
potential due to the Darwin term, 

A t (x) = -<S>(x) + ^V-E(x). (7.10) 

Hence, somewhat unexpected, the Pauli equation also transforms covariantly under 
local SU(2) S transformations 

*(x)' = g(x)*(x), W^x)' = g(x)(W^x) + d,)g{x)l (7.11) 

It should be pointed out that SU(2) S is not a gauge symmetry in the usual sense. 
In particular, the non-Abelian vector potential W^(x) is not an independent degree 
of freedom, but just given in terms of the external electromagnetic field strengths 
E(x) and B(x). The local SU (2) s symmetry is related to the global spatial rotations 
discussed before. In particular, global SU(2) S transformations take the form 

*(x)' = g*(x), W^x)' = gW,(x)g\ (7.12) 

which, for example, implies 

B(xY = O t B(x), (7.13) 

where the resulting 3x3 rotation matrix O G 5*0(3) is again given by eq.(7.5). In 
contrast to a full spatial rotation, a global SU(2) S transformation does not rotate 
the argument x of the magnetic field to Ox. Also the potentials $(x) and A(x) are 
unaffected by the global SU(2) S symmetry. Consequently, the SU(2) S symmetry is 
inconsistent with the relations of eq.(7.2). Despite this, the local SU(2) S symmetry 
of the Pauli equation, which will be inherited by the Hubbard model and by the 
effective theory, dictates how low-frequency external electromagnetic fields are to be 
included in those theories. The high-frequency internal electromagnetic fields (for 
which eq.(7.2) is essential) are integrated out in the effective theory and thus do 
not spoil the symmetry. The local SU(2) S structure implies that in non-relativistic 
systems spin plays the role of an internal quantum number analogous to flavor in 
particle physics. 
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7.2 The Hubbard Model in an External Electromagnetic 
Field 



In the next step we want to couple external electromagnetic fields to the Hubbard 
model. The Frohlich-Studer SU(2) S symmetry of the Pauli equation determines how 
to do this. One must simply use SU(2) S ® U{1)q covariant derivatives with E(x) 
and B(x) playing the role of non-Abelian vector potentials for SU(2) S . Since the 
Hubbard model is defined on a spatial lattice, it is natural to construct corresponding 
SU(2) S <S> U(1)q parallel transporters U x ^ connecting neighboring lattice sites x and 
x + i, 

U X i — VeKp[ ds WAx + si)] explie ds AAx + si)]. (7.14) 
Jo Jo 

Here V denotes path ordering along the link. Under local SU(2) S transformations 

the parallel transporter transforms as 

U' Xti = g(x)U x>i g(x + iy, (7.15) 

while under U(1)q gauge transformations one has 

Q U X;i = exp(iuj(x)) U x , t exp(-iu(x + i)). (7.16) 

The Hubbard model Hamiltonian coupled to external electromagnetic fields then 
reads 

h\u\ = -tj2(4u x ,c x+i + <i./r,,,r,) + \ E(4cx - 1) 2 - ^E(4^ - 1), (7.17) 

x,i x x 

and the corresponding Schrodinger equation takes the form 

iD t V = H[U\V. (7.18) 
Here ^ is the multi-particle wave function and the covariant derivative is given by 

D t = d t + iJ^Wtix) ■ S x + eA t (x)Q x ]. (7.19) 

X 

It should be noted that the Zeeman coupling n e B(x) ■ S x enters the Hubbard model 
through D t , while the spin-orbit coupling appears in the non-Abelian SU(2) S part 
of the parallel transporter U x ^. 

In the Hilbert space of the theory local SU(2) S ® U(1)q transformations are 
implemented by unitary operators 

V = exp(^77(x) • S x ), W = exp(tJ2^(x)Q x ), (7.20) 

X X 

such that 

c' x = V ] c x V = exp(iff(x) ■ ^)c x = g{x)c x , g(x) G SU(2) a , 

Q Cx = W^c x W = exp(iuj(x))c x , exp(iu(x)) e U(1) Q . (7.21) 
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Together with eqs.(7.15) and (7.16) this implies that under the local transformations 
the Hamiltonian transforms as 

H[W] = VH[U]V\ H[ Q U] = WH[U]W\ (7.22) 

Similarly, one obtains 

D' t = VD t V\ Q D t = WD t W\ (7.23) 
such that the Schrodinger equation indeed transforms covariantly when one uses 

m> = V^, Q ^ = W^. (7.24) 

7.3 External Electromagnetic Fields in the Effective Theory 
for Magnons and Charge Carriers 

The couplings of magnons to external electromagnetic fields have been investigated 
in detail in [79]. Again, the Frohlich-Studer symmetry is crucial and one obtains 

S[e, Wg = J d 2 x dt ^ (D£ ■ D t e+ ^D t e- D t ij , (7.25) 

with the covariant derivative 

D^{x) = <9 M e(x) + e(x) x W^x). (7.26) 

Since magnons are electrically neutral, one may expect that they do not couple 
directly to the electromagnetic vector potential A^(x). Still, as discussed in [79] the 
issue is potentially non-trivial because there is a Goldstone-Wilczek current 

tf W ( x ) = e» • [DAx) x D p e(x) + W vp (x)], (7.27) 

with the non-Abelian field strength given by 

W^x) = dpW^x) - ajfrpix) - W^x) x W v (x). (7.28) 

The Goldstone-Wilczek current is an SU(2) S gauge-invariant extension of the baby- 
Skyrmion current of eq.(3.45) and is also topologically conserved, i.e. d^j^ w = 0. 
Hence, one may be tempted to add a Goldstone-Wilczek term j^ w (x)A fl (x) to the 
Lagrangian. However, just like the Hopf term, the Goldstone-Wilczek term breaks 
R, T, and T' and is thus forbidden in the present case. 

Using the P(x) notation, in the presence of external electromagnetic fields the 
action of eq.(7.25) is given by 

S[P, Wg = J d 2 x dt Ps (Tr[DiPDiP] + ^Tr[D t PD t P]j , (7.29) 
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where the SU(2) S covariant derivative is denoted by 

D^Pix) = d,P(x) + [W„(x),P(x)}. (7.30) 

As a consequence of the Frohlich-Studer symmetry, the action of eq.(7.29) is invari- 
ant even under local SU(2) S transformations 

P(x)' = g{x)P{x)g{x)\ W^x)' = g(x)(W^x) + d^)g{x)K (7.31) 

Let us now discuss how external electromagnetic fields enter the fermionic part 
of the effective action. As a rule, ordinary derivatives must be replaced by covariant 
ones. This is the case also in the construction of the composite vector field which 
now takes the form 

v^x) = u{x)D tx u{x^ = u(x)[d^ + W^x)]u(x)l (7.32) 

Under the local SU(2) S symmetry the field u(x) transforms as 

u(x)' = h{x)u(x)g(x)\ (7.33) 

such that 

v^x)' = h(x)u{x)g{x) ] [d ll + g{x)(W^{x) + d^)g{x) ] ]g{x)u(x) ] h{x) ] 
= h(x)u(x)[d^ + W^x)}u(x)^h(xf 

= h(x)(v^x) + d^h(x)l (7.34) 

This is exactly the same transformation behavior as for the global SU(2) S transfor- 
mation of eq.(3.26). In particular, this implies that the U(l) s covariant derivative 
= + ?f^(x)(T3 need not be modified when SU(2) S is turned into a local sym- 
metry. Of course, according to eq.(7.32), v^x) now contains the electromagnetic 
fields E(x) and B(x) through the non-Abelian "gauge" field W^(x). Due to the 
local U(1)q symmetry, the covariant derivatives still need to be extended to 

D^ B = d^ A > B + ivl{x)a^ A > B + ^ A ' B ieA^x)a 3 , 

D^ A ^ = d^ A ^ - m A ^tvl{x)a, - ieA,{x)a^ A ^. (7.35) 

It should also be noted that the low-energy effective theory is not necessarily just 
minimally coupled. In particular, the field strengths F^ v (x) = d^A^x) — d^A^x) 
and W^ u (x) may also directly enter the low-energy effective theory. 



8 Conclusions 

We have constructed a systematic low-energy effective field theory describing the 
interactions of magnons with charge carriers doped into an antiferromagnet. A key 
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ingredient for constructing the effective theory are symmetry considerations. The 
effective theory makes model-independent predictions for magnon-magnon, magnon- 
hole, and magnon-electron scattering. It also determines the long-range magnon- 
mediated forces between electrons or holes. Although these would be highly non- 
trivial non-perturbative issues from the point of view of Hubbard-type models, in 
the framework of the effective theory they can be understood quantitatively by per- 
turbative analytic calculations. More ambitious non-perturbative questions might 
also be within reach of the effective theory. Such questions include the quantitative 
understanding of the Mott insulator state, the reduction of the staggered magnetiza- 
tion upon doping, the formation of a spiral phase, or the systematic investigation of 
dynamical mechanisms for the preformation of electron or hole pairs in the antifer- 
romagnetic phase. In particular, magnon exchange — the analog of pion exchange 
in nuclear physics — suggests itself as a relevant mechanism. 

Before one can do loop-calculations in the effective theory, one must establish 
a consistent power-counting scheme. This has originally been done for pion chiral 
perturbation theory [42], and carries over to magnon chiral perturbation theory in 
a straightforward manner. When charge carriers are included, the issue must be 
reconsidered. The same was true for baryon chiral perturbation theory of pions 
and nucleons. In the baryon number B = 1 sector a consistent power-counting 
scheme enabling a systematic loop-expansion of the effective theory was established 
by Becher and Leutwyler [73]. It is to be expected that this scheme can be ex- 
tended to the low-energy theory of magnons and charge carriers developed here. 
The systematic power-counting in sectors with B > 2 still is a controversial issue in 
baryon chiral perturbation theory. The Weinberg power-counting scheme [52] seems 
to work in most (but not necessarily in all) cases. Its relation to the alternative 
Kaplan-Savage- Wise scheme [53] should be clarified further [59,61]. In light of the 
experience with effective theories for the strong interactions, one should hence ex- 
pect the issue of power-counting to be non-trivial in sectors with two or more charge 
carriers. 

Even when the extra SU(2)q symmetry is imposed, in the fermion sector the 
effective theory has a large number of low-energy parameters. There are two rest 
mass parameters Mi and M 2 as well as two kinetic mass parameters M[ and M' 2 for 
the fermions, four coupling constants A 1; A 2 , K ± , and K 2 for fermion-one- magnon 
vertices, two coupling constants Aq and N 2 for fermion-two-magnon vertices, five 
4-fermion coupling constants G\, G 2 , G§, two 6-fermion couplings Hi and H 2 , 
and finally one 8-fermion coupling /. If only the U(1)q symmetry is imposed there 
are even more parameters. The large number of a priori undetermined low-energy 
parameters is the price one has to pay for the universality and model-independence 
of the effective theory. Only in this way the low-energy physics of any arbitrary 
cuprate antiferro magnet can be captured by the effective theory. Of course, due to 
the rather large number of parameters, the predictive power of the effective theory 
is somewhat limited. Still, only a few parameters enter in some relevant physical 
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quantities. For example, the one-magnon exchange potential between charge carriers 
depends only on certain combinations of the fermion-magnon couplings A 1; A 2 , K ± , 
and K 2 . Also, for example, the details of the short-range 4-, 6-, and 8-fermion 
couplings are not expected to be essential for identifying potential mechanisms for 
preforming electron or hole pairs in the antiferromagnetic phase. It is interesting 
to note that the low-energy effective theory of the t-J model, in which electrons 
are excluded beyond half-filling and holes are the only charge carriers, has a much 
smaller number of low-energy parameters. In that case, there are only one rest 
mass parameter M, one kinetic mass parameter M', two coupling constants A and 
K for hole-one-magnon vertices, one coupling constant A" for a hole-two-magnon 
vertex, and one 4-fermion coupling constant G. It would be interesting to perform 
numerical simulations of the Hubbard or t-J model in order to determine the values 
of the corresponding low-energy parameters by comparison with calculations in the 
effective theory. For example, in the t-J model one can determine the parameters M, 
M', A, K, and A" from simulations in the one-hole sector, while the determination 
of G requires computations in the two-hole sector of the Hilbert space. 

It should be pointed out that, as it stands, the effective theory is applicable only 
at small doping, i.e. for small /i. This is sufficient for understanding the long-range 
forces between electrons or holes in the antiferromagnetic phase. It should also allow 
a quantitative investigation of the reduction of the staggered magnetization upon 
doping. However, in order to enter the high-temperature superconducting phase 
itself, if this is at all possible within the effective field theory presented here, larger 
values of ji will be necessary. Once \i becomes large, it sets a new scale which 
must be taken into account in the power-counting. However, most important, the 
symmetry considerations of the present paper still apply in that case as well. 

Some of the most interesting questions one can address in the framework of the 
effective theory may require non-perturbative calculations. While in some cases 
such calculations can be performed in the continuum, in others they may require a 
non-perturbative regularization of the effective theory. In [92] the effective theory 
of pions and nucleons was regularized on a space-time lattice in order to address 
non-perturbative questions concerning the strong interactions. It may also be useful 
to formulate the effective theory of magnons and charge carriers on the lattice. For 
example, it would be interesting to investigate if the effective theory is more easily 
solvable by numerical simulation than the standard Hubbard-type models. 

We like to emphasize again that effective field theory also allows us to include 
phonons in addition to magnons. This may shed light on more complicated potential 
mechanisms for Cooper pair preformation which involve both magnon and phonon 
exchange. It is interesting to construct such an effective theory. In particular, the 
Galilean (or even Poincare) symmetry is then non-linearly realized. 

To summarize, low-energy effective field theory is a powerful tool that has several 
advantages compared to the direct use of microscopic models. First, it is model- 
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independent and provides universal predictions. Material-specific details of the un- 
derlying microscopic system enter the effective theory only through low-energy pa- 
rameters whose values can be determined by comparison with experiments or with 
numerical simulations. Second, and most important, the electrons or holes of the 
effective theory are quasi-particles whose long-range forces are weak and calculable 
in perturbation theory. This is a significant advantage compared to calculations 
in microscopic models of strongly correlated electrons which are necessarily non- 
perturbative. While it is practically impossible to reliably determine the long-range 
forces between charge carriers from Hubbard-type models, in the effective theory 
the calculation of the one-magnon exchange forces is straightforward and presently 
in progress. It is very interesting to ask if these forces will provide a potential 
mechanism for the preformation of electron or hole pairs. In any case, we propose 
the systematic low-energy effective field theory approach as a better compromise 
between calculability and predictive power than the one offered by Hubbard-type 
models. Effective field theory sheds new light on the dynamics of charge carriers 
in antiferromagnets, and there is hope that it may even be applicable to the high- 
temperature superconductors themselves. 
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A Electron-Hole Representation of the Hubbard 
Model Operators 

For U ^> \t\ the Hubbard model at half-filling reduces to the antiferromagnetic 
quantum Heisenberg model. In contrast to the Heisenberg ferromagnet, the ground 
state of the antiferromagnet is not known analytically. In particular, the naive Neel 
state 



with all spins down on the even sublattice A and all spins up on the odd sublattice 
B is not an eigenstate of the Hubbard Hamiltonian. Still, we use this state in order 
to define electron and hole operators. For even sites we then find 



w = n 4 n 4 T io), 




xeA xeB 



c x] \N) = 0, 4JA0 = 0, xeA. 




46 



Correspondingly, creates an electron, while c x ± creates a hole. Hence, just like a 
relativistic Dirac spinor, the SU(2) S spinor 

-{%)-{%)> ** A ' (A ' 3) 

consists of a particle annihilation operator a x ^ in the upper component and a hole 
creation operator b x + in the lower component. Note that the annihilation of an 
electron with spin down via c x ± corresponds to the creation of a hole with spin up 
via b x y Similarly, on the odd sites one has 

c xl \N) = 0, 4 T |iV) = 0, xeB. (A.4) 

In this case, creates a particle, while c x ^ creates a hole and we write 

*= (Z ) = (£)• < A - 5 > 



B Removal of Non- Canonical Terms by a Field 
Redefinition 



The most general «S?7(2)Q-breaking but ?7(l)g-symmetric terms containing one co- 
variant time-derivative are given by 

-Tr[* A *D t * A + ^D t V B ] + -Tr[^ At a 3 D t V A a 3 - ^ a 3 D t V B a 3 ] 
+ -Tr[^ At A^ B ^s + ^ B] D t m A a 3 ] 

=<^>(2g) + <^>(S- (Bl) 

Here ip A ' B {x) results from a field redefinition that diagonalizes the matrices in the 
previous expression. Only the term proportional to a contains the standard form 
^ A ^d t ip A + ip A ^d t ip A + ip^d t ip B + ip B ^dfip 3 . which implies canonical anticommutation 
relations between fermionic creation and annihilation operators in the Hamiltonian 
formulation. The non-canonical terms (proportional to b and c) can be removed by 
an appropriate field redefinition 
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A± = a ± VWTd>, U ± =(^* v ^ ) . (B.2) 



ty X 

Here U± are unitary matrices with 1,7 gR which obey 

a ± b c \ f / A-t 



It is straightforward to show that the redefined fields ip± B (x) have the same 
symmetry properties of eqs.(5.4) and (5.6) as the original fields tp± B (x). Under the 
SU(2) S symmetry the original fields transform as 

4>£' B (x)' = exp{±ia{x))fe B {x), (B.4) 

and after the field redefinition again 

fe(x)' = yfx^[Xfe(x)' ±Yi> B {x)'\ 

= exp(±ia(x))J\2[Xfe(x) ± Yip B (x)\ = exp(±ia(x))fe(x) , 

fetx)' = ^[^Yfe(xY + X^ B {x)'\ 

= exp(±ia(x))^/\^[TYfe(x) + X$g(x)\ = exp(±ia(x))$ B (x). (B.5) 

Similarly, under the U(1)q symmetry the original fields transform as 

Q fe B (x) = exp(iu)fe B (x), (B.6) 

and again 

Q fe(x) = JH[X Q fe(x)±Y Q 4>l{x)} 

= exp(iuj)^J\2[Xfe(x) ± Yip B (x)] = exp(iu)$±(x), 

mix) = ^[TY Q fe(x) +X q ^j b ( x )} 

= exp(iu)^[TYfe(x) +X^ B (x)} = exp(iuj)^ B (x). (B.7) 

Under the modified displacement symmetry D' one has 

D 'fe B (x) = ±^ B ' A (x), (B.8) 

and after the field redefinition one again obtains 

D 'fe(x) = y/\^[X D 'i)i (x) ± Y D '^ B (x)} 

= ±^[X^j b (x) ± Y^(x)} = ±i> B (x), 



D 'i, B {x) = ^\ T [ T Y D 'fe(x)+X D '^ B (x)] 

= ±J\^[tY^ b (x) + X^(x)} = ±$$(x). (B.9) 
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Since the displacement symmetry D is a combination of D' and SU(2) S it also 
maintains its original form. The same is true for the discrete symmetries O and R. 
Finally, under the modified time-reversal V the original fields transform as 

T '^t' B (x) = —ip± B ^(Tx), (B.10) 

such that 

= -yf\^[Xrj}^{Tx) ± Y^(Tx)\ = -^(Tx), 
T 'mx) = ^[TY T '^(x)+X T '^(x)} 

= ~\[K[TY'il)i\Tx)+X^l\Tx)\ = -$ B \Tx). (B.ll) 

As a combination of T' and SU(2) S the time-reversal symmetry T also maintains its 
original form after the field redefinition. The only symmetry that does not maintain 
its original form is SU(2)q. This is no problem since the non-canonical terms can 
arise only when the SU(2)q symmetry is explicitly broken down to U(1)q and is 
hence no longer a symmetry of the theory. 

Since the redefined fields transform exactly like the original ones, the terms in 
the effective Lagrangian take exactly the same form as before. Hence, it is indeed 
justified not to include the non-canonical terms in the effective Lagrangian. 
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